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CHAPTER  I 

INTRODUCTION 

This  paper  deals  with  the  microwave  induced 
heating  in  thin  multiple  conductive  coatings.  An  under¬ 
standing  of  this  phenomenon  is  needed  so  an  optimum 
coating  may  be  designed  for  the  thermographic  detection 
of  surface  currents.  This  process  uses  a  thermographic 
camera  to  detect  temperature  variations  on  an  object's 

surface.  These  temperature  variations  are  caused  by 
,  2 

joule  heating  (I  R)  due  to  surface  currents.  Thus,  it 
is  possible  to  relate  surface  currents  to  surface 
temperature  variations.'  Once  this  is  done  the  thermo¬ 
graphic  system  can  measure  surface  currents  on  an  entire 
object  very  quickly  and  relatively  easily.  However, 
there  is  still  a  significant  limitation  to  this  type  of 
measurement  system.  The  limitation  is  that  we  measure 
a  single  scalar  quantity  (temperature)  at  each  point  on 
the  surface;  hence,  we  only  obtain  the  magnitude  of  the 
total  surface  current  at  that  point  and  have  no  phase  or 
direction  information. 

Nonetheless,  current  aunplitude  Information  alone 
is  a  valuable  piece  of  information.  For  exaunple,  in  the 


2 


study  of  Electromagnetic  Pulse  (EMP)  phenomena  we  might 
want  to  know  where  the  greatest  current  density  ampli¬ 
tudes  are  located  on  an  aircraft  to  aid  in  the  proper 
placement  of  cable  bundles.  Or  possibly  in  radar  cross 
section  studies,  we  might  want  to  know  where  the  greatest 
current  densities  exist  before  a  design  or  numerical 
solution  is  attempted.  And  lastly,  measured  surface 
current  amplitude  information  would  be  extremely 
valuable  in  the  verification  of  numerous  numerical 
techniques  designed  to  provide  amplitude,  phase,  and 
direction  information.  This  thesis  which  relates  coating 
temperature  to  incident  microwave  power  levels  is 
organized  in  the  following  manner. 

The  second  chapter  begins  with  a  discussion  of 
electromagnetics  in  general.  It  considers  the  phasor 
form  of  Maxwell's  equations  in  a  conductive  medium  and 
then  reviews  their  simultaneous  solution.  The  result 
yields  the  type  form  the  wave  must  take  as  it  propagates 
in  the  medium  and  discusses  some  of  its  characteristics. 
Particular  areas  of  discussion  include  attenuation  in  a 
medium  and  the  complex  form  of  Snell's  law. 

The  third  chapter  deals  exclusively  with  the 
solution  of  a  single  boundary  electromagnetic  problem. 

The  boundary  is  a  semi- inf inite  plane  located  at  the 
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zaO  interface  between  air  and  a  conductive  medium.  The 
solution  is  completed  for  various  incident  angles  and 
for  polarizations  both  parallel  and  perpendicular  to  the 
plane  of  incidence.  The  single  boundary  problem  provides 
two  useful  Insights  to  our  overall  problem.  First  »  it 
is  done  in  a  general  manner  so  the  seune  overall  approach 
is  used  in  the  solution  to  the  N-layer  system  later  in 
chapter  four.  Second,  the  results  may  be  used  to  provide 
information  concerning  the  theoretical  infrared  emissiv- 
ity  of  the  coatings  that  are  to  be  placed  on  the  various 
object  surfaces.  The  rate  and  direction  of  Infrared 
emission  are  both  directly  proportional  to  this  emissiv- 
ity  factor.  Therefore,  our  detection  effectiveness  is 
very  dependent  on  surface  emissivity.  Various  plots  are 
presented  which  illustrate  the  reflectivity  and  emissivity 

7 

for  materials  with  conductivities  of  0  to  5x10  mhos/m. 

Chapter  four  presents  a  detailed  solution  for 
the  N-layer  electromagnetic  interaction  problem.  The 
problem  considers  both  parallel  and  perpendicular 
polarizations  at  a  particular  incident  angle.  Each  layer 
is  considered  to  have  a  particular  electrical  conductiv¬ 
ity.  The  solution  is  developed  in  terms  of  an  expanded 
matrix  in  a  manner  similar  to  that  used  in  chapter 
three.  A  computer  program  is  included  in  appendix  D 
that  will  do  the  solution  for  up  to  a  10  layer  system 
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and  display  the  results  with  a  Hewlett-Packard  9845B 
mini -computer.  Included  is  an  exaunple  of  the  plot 
obtained  for  a  three  layer  system  with  the  layers  having 
conductivities  of  .5,  1.0,  and  1.5  mhos/m  respectively. 

Of  note  in  this  development  is  that  there  are  no 
approximations  used  other  than  assuming  semi-infinite 
sled)  dimensions. 

In  the  next  chapter  we  tie  the  electromagnetic 
interaction  to  the  thermodynamics  of  the  problem.  Thus, 
for  a  given  layer  configuration  we  relate  surface  temper¬ 
ature  to  a  given  microwave  input  power.  Initially  in 
the  development,  a  simplified  model  is  presented  that 
allows  us  to  estimate  the  transient  response  of  the 
system.  Then  somewhat  later  the  steady  state  thermal 
characteristics  of  the  system  are  analyzed.  In  this 
analysis  we  assume  there  are  no  heat  losses  from  other 
than  the  z«0  interface;  that  is,  we  assume  one  of  the 
layers  below  the  conductive  coating  to  be  a  perfect 
insulator.  At  the  z»0  interface  we  consider  both  infra¬ 
red  and  convective  heat  losses.  The  steady  state 
solution  to  the  heat  equation  yields  a  transcendental 
equation  that  must  be  solved  numerically.  A  computer 
program  is  written  that  does  this  solution  and  plots 
surface  temperature  versus  electrical  conductivity  and 
coating  thickness  on  a  three  dimensional  plot.  Plots 
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are  illustrated  for  typical  two  and  three  layer  systems. 

As  in  any  theoretical  development »  we  must  have 
experimental  verification  before  the  results  may  be  used 
with  any  confidence.  Chapter  six  provides  experimental 
verification  of  the  computer  results  for  various  layer 
configurations.  The  three  primary  sources  of  error  in 
the  experiments  resulted  from  the  lack  of  semi -infinite 
planes  to  irradiate  with  microwaves;  thus»  samples  much 
smaller  than  a  wavelength  were  used  to  avoid  resonance 
problems.  Also,  the  problem  of  assuming  zero  heat 
conduction  losses  from  the  interface  was  partially 

solved  by  placing  the  samples  on  a  styrofoam  substrate. 
Finally,  the  third  major  source  of  error  resulted  from 
the  convective  heat  transfer  coefficient  used  which 
itself  was  an  empirical  value.  See  for  example  Holman's 
book.  Heat  Transfer.  Nonetheless,  predicted  surface 
temperatures  were  very  close  to  the  measured  values. 

Chapter  seven  is  a  mostly  qualitative 
discussion  of  the  relationship  of  our  one -dimensional 
model  to  the  more  complicated  two  or  even  three- 
dimensional  "real  world"  experiments.  In  particular, 
the  two  areas  discussed  involve  the  validity  of  using 
a  small  finite  shape  in  the  experimental  verification 
of  Chapter  VI  and  the  question  of  "nearest  neighbor" 
significance  in  the  actual  measurement  of  surface 
currents  on  electrically  large  targets.  Theoretical 
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as  well  as  experimental  results  are  presented  that  help 
establish  the  validity  and  usefulness  of  our  one- 
dimensional  model  as  applied  toward  the  '*real  world" 
situation. 

In  Chapter  VIII  recommendations  and  conclusions 
are  made  in  regards  to  what  to  use  as  an  optimiim  coating 
in  the  thermographic  detection  of  microwave  induced 
surface  currents.  Basically,  there  are  two  situations 
that  must  be  considered.  The  first  considers  geometry 
only.  If  we  are  interested  only  in  the  geometrical 
aspects  of  a  problem  such  as  studying  the  currents  on  an 
entire  aircraft,  the  simpler  approach  is  to  build  a 
model  of  foam  and  then  coat  it  with  a  thick  (two  skin 
depths)  layer  of  material  with  a  conductivity  of  from 
300  to  750  rahos/m.  Hence,  the  currents  we  measure  are 
the  currents  in  the  coating  itself;  we  simulate  a 
highly  conductive  object  with  one  of  lesser  conductivity. 
The  other,  and  more  difficult,  situation  is  one  in 
which  the  material  an  object  is  made  matters  and  cannot 
be  approximated.  This  situation  would  occur  if  several 
different  materials  were  used  in  the  construction  of  a 
particular  shape  (For  example,  there  might  be 
ferromagnetic  materials,  copper,  aluminum,  and  composites 
all  located  on  the  underside  of  an  aircraft  of  space¬ 
craft.).  Here  we  would  like  to  determine  how  these 
different  materials  behave  together.  Therefore,  we 
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place  our  coating  directly  on  the  object  itself.  By 
appropriate  design  of  the  coating,  we  can  detect  the 
current  distribution  under  our  coating. 

Various  coating  options  are  described  in  the 
final  chapter  as  well  as  their  advantages  and  disadvan¬ 
tages.  The  various  appendices  contain  the  experimental 
arrangements  for  measuring  such  material  characteristics 
as  electrical  conductivity  and  permittivity.  Also 
included  are  all  the  major  computer  programs  used  in  any 
of  the  numerical  analysis  with  the  9845B  computer.  This 
computer  uses  a  Hewlett-Packard  enhanced  BASIC  language 
that  would  be  easily  modified  for  use  on  a  computer 
using  either  Fortran  or  Pascal. 

Overall  then,  this  paper  begins  by  presenting 
the  development  of  a  one-dimensional  model  for  the 
electromagnetic  interaction  with  a  system  of  N  layers 
with  differing  electrical  characteristics.  We  then 
use  this  model  in  a  thermodynamic  analysis  to  arrive 
at  the  steady  state  surface  temperature  of  the  system 
of  layers.  Thus,  we  arrive  at  a  model  which  predicts 
the  equilibrium  surface  temperature  of  a  system  of  N 
discrete  layers  in  the  presence  of  electromagnetic 
radiation.  This  model  is  verified  experimentally  and 
its  applicability  to  two-dimensional  targets  is 
discussed.  Finally,  recommendations  are  made  in  regards 
to  particular  coating  design. 
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CHAPTER  II 

GENERAL  ELECTROMAGNETIC  ANALYSIS 


In  electromagnetic  analysis  we  will  study  the 
simple  one  boundary  problem  and  then  progress  to  multiple 
boundaries.  Throughout  this  analysis  we  will  assume  a 
time  dependence  for  E  and  H  of  the  form  We  will 

first  calculate  the  form  the  field  must  take  by  looking 
at  the  solution  of  Maxwell's  equations  in  a  conductive 
medium. 


In  a  charge  free  conductor  and  with  the  assumed 
time  dependence,  Maxwell's  equations  may  be  written  as 


follows: 

7  X  E  =  jujvi  H  (1) 

V  X  H  =  (o-  jwe)  E  (2) 

V  •  E  *  0  (3) 

V  .  H  =  0  (4) 

Solving  equations  (1)  through  (4)  simultaneously^^  we 
can  arrive  at  the  vector  wave  equation 

V'  E  +  E  *  0  (5) 


where 

is  defined  by 
0  0 

=(i)^UE  +  j  uyo 


(6) 
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(14) 
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From  this  point  on  we  will  omit  the  e~  time 
dependence,  since  it  will  occur  in  all  of  our  terras,  and 
proceed  with  the  more  familiar  phasor  notation.  Thus  we 


will  write  for  the  electric  field 


E 


(15) 


This  is  a  shorthand  notation  and  it  is  to  be  understood 
that  the  is  always  present  even  though  it  is  not 

written , 

Lastly,  it  is  advantageous  to  develop  the  complex 
Snell's  law.  Consider  the  single  interface  below 
separating  two  conductive  materials. 


EM 

Wave 


Figure  1:  Plane  Wave  Interaction  with  a  Single  Interface 
We  know  the  wave  phase  must  be  the  same  on  each  side  of 
the  boundary  at  any  given  point  on  the  boundary;  other¬ 
wise,  the  wave  would  "tear"  and  dispersion  would  occur. 

1 2 

This  phase  requirement  may  be  satisfied  by  having 
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Again  referring  to  the  drawing,  we  have 

e2(cos0  X  +  sine  y)T  s  ^2  x  +  sini|i  y)  t  (17) 

Since  r  =  y  y  defines  the  interface,  the  result  is 

5-^  sine  =  62  sin<i>  (18) 

Thus,  sin4>  is  in  general  complex  and  is  given  by 


+  jYi  . 

sin^  =  ^  sine  =  sine 


(19) 


We  may  also  write  the  expression  for  cosij'  by  recogniz- 

2  2  13 

ing  that  for  any  complex  z,  sin  z  +  cos  z  =  1  ;  thus, 


COSdi  = 


2  2 

ai  -  Yi  +  j  201  y-,  , 

1 - = - k - sin'^e 

2  2  ^  , 

«!  -  Y2  *  3  ^^2  ^2 


(20) 


We  can  loolc  at  Snell's  Law  in  the  case  of  an 
air/good  conductor  interface.  If  the  rtiedium  on  the 

Ott 

left  is  air,  reduces  to  —  ,  where  >■  is  the  free 

space  wavelength.  In  the  good  conductor  we  have 

1  2Jl  _ _ 

a=  Y=  —  ;  then  we  have  sin(Ji  =  (-t — )(a-  +  jY^  )sine  (21) 

6  i  L 

with  algebraic  manipulation  this  reduces  to 


sincfi 


1 

2c/y 2^2 


(1  -  j)  sine 


(22) 


In  the  good  conductor  limit  we  have  therefore, 

sin<ti  =  0  (23) 

which  implies  <ti  =  0 

Hence,  regardless  of  the  incident  angle  e,  the  wave  will 
propagate  approximately  normal  to  the  surface  after 
entering  medium  2. 


12 


CHAPTER  III 

SINGLE  BOUNDARY  ANALYSIS 

With  the  basic  equations  developed,  we  can  now 
do  the  actual  field  calculation.  We  will  begin  with  the 
single  boundary  which  will  provide  two  important  pieces 
of  information.  First,  it  will  provide  needed  insight 
for  the  solution  of  the  more  complex  multi-layer  problem, 
and  second,  it  will  provide  a  method  of  calculating  theo¬ 
retical  infrared  spectral  emissivities.  This  information 
will  be  very  valuable  in  the  thermodynamic  analysis  of 
the  multi-layer  problem  later  in  the  development. 

We  begin  by  considering  the  diagram  below  in 
which  we  have  an  incident  electromagnetic  wave  from  the 
left.  Parallel  incidence  is  illustrated. 


Figure  2:  Plane  Wave  Interface  on  a  Single  Boundary 
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From  the  previous  development  we  know  that 


a  - 

(24) 

y  -  7^/-! 

(25) 

62““  +  vIY  ^0=^ 

(26) 

sin*  *  B  Sin8 

(27) 

,  =  /i  o2  sin'^e  j  2aY) 

V  °  r  2.  2.2 

V  (a  +Y  J 

(28) 

We  will  need  the  real  and  imaginary  parts  of  co8i<> 
later  so  it  is  advantageous  to  calculate  them  now.  We 
begin  by  letting  cos<fi  =  x  +  j  y  where  x  and  y  are  de¬ 
fined  to  be  real .  Thus , 

cos^'P  =  -  y^  +  j  2xy  ^29) 

Equating  real  and  imaginary  parts ,  we  obtain 
two  simultaneous  equations; 


x2  -  .  1  - 

(o^  + 

*  fi2  sin^e  c»v 

xy  -  - 7 — 7—7 


(30) 


(31) 


Solving  (30)  and  (31)  by  substitution,  yields  the 
following 

X  -  \/p+  \/p^  + 


(32) 


y?  +  y 


p2  +  q2 


y 


(33) 
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where 


sin^e (g^-Y^) 

2(g2  +  Y^)^ 


sin  9  g  Y 

(g2  .  y2)2 


(34) 

(35) 


In  the  good  conductor  limit  where  (— )>>  1  these  reduce 

^  toe 

to 


For  example  if  the  incident  frequency  is  2.5  GHz  and 
a  =  10  mhos/m  we  have  cos<t>  =  1.00002  +  j(.007)  which 
indicates  that  <(>  may  be  considered  to  be  approximately 
zero  for  values  of  o>  10.  In  some  numerical  calculations 


this  could  vastly  shorten  the  computer  time. 

By  considering  the  geometry  of  the  problem,  we 
may  write  the  wave  vectors  as  follows; 

=  6o  (sine  X  cose  z)  (38) 

^2  =  6o  (sine  X  -  cose  z)  (39) 

^  A  ^ 

tSj  =  62  (sin'll  X  +  cos<)i  z)  (40) 

Since  we  are  considering  an  infinite  plane,  we  are 
interested  only  in  the  z  dependence  of  the  fields;  thus 
we  may  ta)ce  r  =  z  z.  We  have 

r  =  9^  z  cose  (41) 

62*  r  =  z  cose 


(42) 
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r  =  (a  +  jy)  z  COS*  (43) 

At  this  point  we  divide  the  problem  into  two  parts  since 
we  must  consider  parallel  and  perpendicular  polarization 
of  the  E  vector  separately. 

E  Parallel 

If  E  is  parallel  to  the  plane  of  incidence  we 
may  write  the  vector  fields  with  reference  to  the  dia¬ 
gram.  We  have 

•  Ej  °°°°  (cose  X  -  sine  z)  (44) 

H  .  ^  E,  ei  y  (45) 

(1)U 

E,  =  E,  e”^  (cose  X  +  sine  z)  (46) 

H,  =  -  !s  E,  e-3  E  co=»  y  ,47, 

^  uiy  ^ 

E3  =  E3  e^  (cos*  X  -  sin*  z)  (48) 

«3  =  e3  ^2  ®  y 

To  satisfy  the  boundary  conditions  at  z  =  0,  we  equate 

the  tangential  field  components  of  E  and  H  there;  thus, 

we  have 

=  E3  cos*  (50) 


Ej^  cose  ♦  E2  cose 

E.  -  ^  E~  “  ^2 

(oy  ^  toy  ^ 


(51) 
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(58) 


As  before,  we  apply  the  boundary  conditions  at  the 
z  =  0  interface  which  yields 
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El  +  E2 


BqCOS0 

ojy 


6nCOS6 

E,  +  -2 -  E- 

•I-  Ojy  ‘ 


COS(ti 

ujyZ 


(59) 

(60) 


Again  these  equations  simplify  and  may  be  put  into 
matrix  form. 


"1  -1 

N> 

- 

El 

X 

s 

^0  cos  9  ^  cos  (J 

^3 

^0  cos  0 

y 

(61) 


Analytic  Solution 

Equations  (52)  and  (61)  may  be  solved  directly 
with  little  difficulty;  however,  since  we  are  develop¬ 
ing  a  general  technique  for  a  multi-layer  system,  it 
is  worthwhile  to  do  the  solution  numerically.  The 
technique  is  relatively  straight  forward  in  that  the 
coefficient  matrices  are  loaded  into  a  computer  along 
with  the  constant  vector.  For  example,  if  we  have  the 
matrix  equation  Ax  =  B,  we  would  load  the  auxiliary 
matrix 


®11  ®ln 

'<  ^22  ! 

> 

•  ^ 

®nl  *nn 


(62) 


into  the  computer.  It  would  be  an  N  x  N  -t-  1  dimensional 
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matrix.  We  then  do  a  Gauas-Jordan  row  reduction  to  put 
the  A  portion  into  identity  format.  The  B  column  will 
now  correspond  to  the  solution  vector  of  x.  This  is 
relatively  easy  to  program  unless  A  happens  to  be  com¬ 
plex.  In  this  case  the  complex  equation  Ax  =  B  may  be 
expanded  into  its  real  and  imaginary  parts  to  yield  a 
matrix  equation  of  only  real  terms  which  may  be  handled 
in  a  straight  forward  manner.  For  example  if  the  set 
of  equations. 


11  *1 

*  ®12  *2  =  ^1 

(63) 

21  *1 

+  ^22  *2  “  ^2 

(64) 

is  complex,  we  may  in  complex  notation  write  a  complex 
matrix  equation  as 


®12 


®21  ®22 


(65) 


Likewise  we  may  write  (63)  and  (64)  in  an  expanded 
format  in  terms  of  its  real  and  imaginary  parts.  If 
we  let  R  and  I  designate  Real  and  Imaginary  respec¬ 
tively,  we  have 

(Raii+jiaii)  (RXi+jlXi)  +  (Raj^2*3laj^2)  <***2*^^^2^  ” 

Rbj^  + 


(66) 
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(Ra2j^+jla2^)  (RXj^+jlXj^)  ♦  (Ra22+jla22)  “ 

Rb2  +  jlb2  (67) 

Equating  real  and  imaginary  parts  in  (66)  and  (67)  we 
have 


RXj^ 

* 

'“ii 

IXi 

+ 

Rai2 

RX2 

— 

1^12 

1x2  =  Rb^ 

(68) 

RXj^ 

+ 

Rau 

IXi 

+ 

Iai2 

RX2 

+ 

^®12 

1X2  =  Ibj^ 

(69) 

RXj^ 

- 

lajj 

IXi 

^*22 

RX2 

- 

’•^22 

1X2  =  Rb2 

(70) 

RXj^ 

+ 

"“21 

IXi 

^®22 

Rx2 

+ 

^®22 

1X2  =  Ib2 

(71) 

or  in  matrix  format  we  have  a  matrix  of  real  terms  only 
written  as  follows 


"“11 

-^“11 

• 

I 

1 

"“12 

"^®12 

»  * 

RXj^ 

»  ^ 
Rbj^ 

'“11 

"“11 

1 

( 

1 

^“12 

"“12 

X 

IX^ 

S 

"“21 

1 

f 

Ra22 

"^^22 

RX2 

Rb2 

"“21 

1 

1 

^^22 

^®22 

1X2 

Thus,  to  numerically  solve  a  2  x  2  complex  matrix  a 
computer  must  reduce  a  4  x  4  matrix,  and  in  general  to 
solve  an  N  X  N  complex  matrix  we  must  generate  a  2N  x 
2N  real  matrix. 

There  are  some  advantages  to  this  approach  in 
that  after  the  reduction  is  complete  we  have  all  our 
variables  written  in  real  and  imaginary  format.  Also, 
the  expanded  matrix  is  easy  to  generate  if  one  parti¬ 
tions  the  original  complex  matrix  and  then  notes  that 
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each  complex  element  goes  into  a  2  x  2  matrix  that  is 
nearly  semetric  except  for  one  sign. 

Wtimeric  Solution 

The  numeric  solution  takes  the  matepial  char¬ 
acteristics  (°2'^2'  ^2^  then  constructs  the  expanded 
auxiliary  matrix,  Para  (I,J)  and  Perp  (I,J)  correspond¬ 
ing  to  the  Parallel  or  Perpendicular  solutions  respec¬ 
tively.  The  program  then  does  a  Gauss-Jordan  elimination 
for  each  of  these  as  it  iterates  through  1°  increments 
from  0  to  90°.  For  each  iteration  a  reflectivity 
coefficient  is  calculated  and  stored.  After  the  routine 
is  complete  the  reflectivity  coefficients  are  plotted 
as  a  function  of  d  for  both  parallel  and  perpendicular 
polarizations  on  the  same  plot.  Additionally  a  plot  of 
spectral  emissivity  (absorbance)  is  calculated.  For 
the  emissivity,  X  is  set  equal  to  3  ym  and  is  calculated 
as  e  =  (1  -  Reflectivity).  The  result  is  then  plotted 
as  e  versus  0.  The  program  written  for  a  Hewlett- 
Packard  9845B  computer  is  listed  in  appendix  C.  Figures 
3  through  10  are  the  microwave  reflectance  and  emissiv¬ 
ity  plots  for  a  dielectric  and  a  conductor  with  conduc- 
tivities  of  0,  3,  500,  and  5x10  mhos/meter  respectively. 


-OOT 


Figure  3;  Reflectance  from  a  Single  Plane  Interface  (o=0) 


Pe  r  p  end  i  c  u  a  n  c  e 


100 


24 


mhos/ 


Parallel  Emissivity 
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Since  we  have  N  separate  layers,  ^2Ti-*-l 
electric  field  of  the  wave  traveling  to  the  right  in 
material  n  and  ^2n+2  wave  traveling  to  the 

left.  As  in  the  single  boundary  case,  we  may  write  the 
field  components  in  each  layer  as  follows; 

E  Parallel 

je  z  COS0 


Ef  =  E^e 


H,  =  Eie-"'o 
1  uy 


O  (COS0Q  X 


sinOjj  z) 


jB  z  COS0 


o  y 


t-  =  E-e'^^o^  (cos0  X  +  sin0  z) 

2  2  O  O 


«2  = 


^o 

ojy 


„  -i6  z  COS0  - 

E^e  o  y 


Ej  =  E^ej^l^  COS0J  (COS0J  x  -  sin0j  z) 


H3  = 


^1 


E  JB^Z  COS01  ^ 


E^  =  E^e  cos9j^  (cosQj^  x  +  sinB^^  z) 

6- 


»4  = 


’1  y.  COS0  y 


a>v) . 


E^e 


(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 


'2n+l 


2n+l 


^2n+l  cose^  (cosBj^x  -  sine^z) 


6 


n 


uy 


n 


‘^2n+l  ®  ^  ^ 


^2n+2  °  ^2n+2 


e'jBnZ  cos0^  ^  +  sin0  z) 

^  n  n 


’2n  +  2 


^n  ^  -iS  z  cos9 

^2n*2  ®  "  y 


c  «  c  16  Z  COS9  ,  .  »  ■  ^ 

^2n*5  ^2n*3  ®  °  °  (cos6^x  -  sin0^  z) 


H 


2n  +  3 


c  cos0^  y 

—  E.^^.e-'  0  o 

mu  2n  +  3 


(81) 

(82) 

(83) 

(84) 

(85) 

(86) 
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Application  of  the  tangential  boundary  conditions 


at  the  interface  yields  the  following: 

E  Parallel 

z  =  0 : 

E,  cos9  +  E_  COS0  =  E,  COS0,  +  E.  cos®-, 
1  o2  o3  14  1 


—  -  E-  —  = 

E,  — 

-  E, 

®1 

y  2  y 

3  yj 

4 

‘^l 

6, z. 

cose , 

«'3  3-1  z. 

COS0 

cose^  e-*  1  1 

1 

+  E4 

cose 

e  1  1 

cose,  e^^2^1 

cose  2 

^  ^6 

cose2 

e-i»2"l 

cose 

(101) 

(102) 


(103) 


E,  !i  -  E,  e-3^^ 


1  = 


E,  ti  “^*2  -  E^  tl  e-j«2"l  =°=*2 
5  ^2  6  ^2 


(104) 


z  = 


r.  ^  z  cose  _  ^  „  -i6  z  cose 

E2n+1  cose^  e-'^^n  n  n  +  E2^^2  e  n  n  n 

E.  cose  ej®o^n  ^°"®o 
Zn+o  o 


(105) 


“n  ei  Vn  ^°=“n  -  E,_,  ^  e'iVn  - 


n 


2n+l 

E  ej^o^n 

‘'2n+3  u  ®  °  " 


'2n+2  n 


n 


n 


(106) 


33 


z  =  z^: 


r,  je  z  cose  i,  -is  z  cose  c  „j6  z  cose^ 

^2n+l  ®  ri  ^  ^2n+2  ®  ^  “  ^2n+3  ®  o  ° 

(111) 

„  ®ncos9„  _jB„z„cos9  ®ncos  -jg  z  cose  _ 

■^211+1  TT-  n  eJ  n  n  n  +  £2^+2  ~  n  e  n  n  n  - 
n  n 

(112) 


0 

ncose 

■^2n+l 

y 

n 

-F 

^ocose 

2n+3 

Thus,  we  have  generated  N  equations  with  N  unknowns  for 
both  parallel  and  perpendicular  incidence.  These  equa¬ 
tions  may  easily  be  consolidated  into  matrix  format;  see 
Figures  12  and  13.  Figure  14  is  an  example  plot  of  a 
three  layer  system  with  layer  conductivities  of  .5,  1,  and 


ignetic  Respons 
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1.5  mho/m  respectively.  See  Appendix  D  for  a  listing  of 
the  computer  program  written  for  a  Hewlett-Packard  9845B 
computer;  it  will  handle  up  to  10  layers. 

It  should  be  pointed  out  that  this  is  not  the 

first  time  the  electromagnetic  N-layer  problem  has  been 

solved.  Hansen  develops  a  procedure  in  which  explicit 

equations  are  derived  for  the  mean-square  electric  fields 

induced  by  plane  electromagnetic  radiation  in  a  single 

boundary,  single  layer,  and  N-layer  system.  His  approach 

is  to  allow  for  a  complex  permittivity  in  the  energy 

absorption  process  rather  than  using  Ohm's  law  explicitly 

■*  in  the  solution  of  Maxwell's  equations.  It  is  for  this 

reason  and  as  a  result  of  the  desire  to  present  a  more 

straight  forward  approach,  that  the  preceding  chapter 

was  developed.  Additionally,  interested  readers  may 

17 

also  refer  to  Wait  for  an  in-depth  discussion  of 
electromagnetic  absorption  by  stratified  media.  Most 
of  Wait's  discussion  is  directed  at  the  electromagnetic 
interaction  occurring  at  the  surface  of  a  stratified 
earth. 
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CHAPTER  V 

THERMAL  ANALYSIS 

The  simplest  analysis  Is  one  that  assumes  a 
steady  state  situation  with  no  heat  losses  at  the  z=z 

n 

Interface.  In  practice  this  Is  not  a  bad  assumption 
since  our  Nth  layer  may  be  made  of  a  low  thermal  conduc¬ 
tivity  material  such  as  one  of  the  many  foams  available. 
Also,  It  has  been  demonstrated  that  for  relatively  thin 
layers  (thickness  less  than  one  millimeter)  steady 
state  Is  reached  very  quickly  (In  most  cases  less  than 
five  minutes  for  10  mW/cm  Incident  power). 

A  rough  estimate  for  the  transient  response  of 
the  system  may  be  made  as  follows.  Consider  the  semi¬ 
infinite  slab  insulated  at  z^L  as  shown  in  Figure  15 


Figure  15:  Three  Layer  Thermal  Transient  Response 

2 

-  total  absorbed  power  (watt/m  ) 

Q  =  total  heat  loss  (Newton  cooling) 

u  *  temperature  of  the  layer  (K)  which  is 
assumed  uniform  throughout 


If  u  is  assumed  constant  throughout  our  thin 
layer  (this  assumes  the  thermal  conductivity,  k,  is 
large  or  the  thickness,  L,  is  small  or  both  for  the 
layer) ,  we  may  write  the  energy  equation  for  the  system 
as  follows: 
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layers  are  heated  at  a  relatively  constant  rate  (same 
u  throughout) ;  hence  the  time  constant  for  the  layered 


system  may  be  approximated  by 


o  . c. L. 

Pill 


(118) 


It  is  important  to  point  out  that  the  above 
development  assumed  a  linear  heat  loss  at  the  z=0  sur¬ 
face.  In  fact,  this  term  consists  of  a  convective  term 
and  a  radiative  (infrared)  term  which  is  in  no  way  lin¬ 
ear  for  large  temperature  increases  (large  P_).  However, 
the  steady  state  surface  temperature  may  still  be  calcu¬ 
lated  using  numerical  techniques  for  this  non-linear  case. 

We  begin  the  steady  state  solution  as  before  by 
writing  the  energy  equation  for  the  system;  thus,  we 


have  absorbed  power,  P  ,  equal  to  convective  losses, 

cl 

Q  ,  plus  infrared  losses,  Q.  .  For  the  two  loss  terms 
con  ir 

^  18,19,20 

we  have 


where 

2 

h  =  convective  heat  transfer  coefficient (watt/m  -K) 
0^*  steady  state  surface  temperature  (K ) 
ambient  air  temperature  (K) 
surface  emissivity 

y  =  Stefan-Boltzman  constant(5.67xl0  ®watt/m^-K^) 
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It  has  been  shown  that  h  for  a  vertical  flat  plate  may 

21 

be  approximated  by 

.  1.42  (0^  -  0  . 

h  =  o  air 


where  H  is  the  plate  height  and  5  is  a  number  that 
falls  in  the  range  .1  <5 <.6.  The  number,  6 ,  is  a 
function  of  such  things  as  atmospheric  pressure,  humid¬ 
ity,  etc.  Furthermore,  P  is  simply  the  incident  power 
minus  the  reflected  and  transmitted  power;  therefore, 
we  may  write. 


It  is  important  to  note  that  it  is  necessaxry  to  consider 
both  convection  and  infrared  radiation  in  the  tempera¬ 
ture  range  near  20C  as  indicated  in  Figures  16  through  18 . 

This  figure  illustrates  the  ratio  of  Q.  /Q  for  a 

ir  con 

surface  temperature  increase  of  IOC  with  an  ambient  air 
temperature  of  20C.  This  ratio  is  nominally  between 
1  and  2  over  this  range  indicating  that  both  processes 
are  significant  in  the  energy  loss  transport  and 
therefore  each  must  be  considered. 


RRTIO  OF  INFRARED  TO  CONVECTION  LOSSES 
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(Uo-Uair)  IN  DEGREES  CENTIGRADE 


Figure  16: 


Infrared  Versus  Convection  Comparison 


RRTIO  OF  INFRRRED  TO  CONVECTION  LOSSES 


(Uo-Uair)  IN  DEGREES  CENTIGRRDE 


Figure  17:  Infrared  Versus  Convection  Comparison 


RRTIO  OF  INFRRRED  TO  CONVECTION  LOSSES 
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5  PLATE  HEl GHT=  . J0  METERS 

I  CONVECTION  EXPONENT=.600 

1  fllR  TEMPERfiTURE  (0=20,00 


Fe=  2 


0 


10 


(Uo— Uair)  IN  DEGREES  CENTIGRRDE 


Figure 


18:  Infrared  Versus  Convection  Comparison 
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Even  though  equation  119  is  relatively  simple 
in  appearance,  it  is  important  to  recognize  that  is 
in  fact  a  function  of  several  variables.  These 
variables  include  the  layer  thicknesses,  relative 
permittivities,  relative  permeabilities,  electrical 
conductivities,  surface  emissivity,  convection 
exponent,  convection  coefficient,  sample  height, 
incident  microwave  power  level,  microwave  frequency, 
and  microwave  incident  angle.  Thus,  in  general  is 
a  function  of  4N+7  variables  where  N  is  the  number  of 
layers  considered.  Even  for  the  case  of  N=l,  we  see 
that  must  be  considered  in  11  space  in  order  to 
analyze  all  the  variables  and  their  interaction  at 
once.  Clearly,  we  can  visually  represent  3  space  but 
not  11;  therefore,  in  order  to  graphically  display  the 
interactions  of  the  more  dynamic  variables- -the  ones 
over  which  we  have  direct  control --we  fix  the  values  of 
such  things  as  incident  angle,  power  level,  convection 
exponent,  etc.,  and  allow  only  coating  thickness  and 
coating  electrical  conductivity  to  vary  which  then 
yield  our  equilibrium  surface  temperature,  U^. 
Topologically  then,  by  allowing  only  two  of  the 
variables  to  vary  we  are  in  effect  looking  only  at  a 
particular  plane  in  4N+7  space.  As  was  indicated 
earlier,  this  is  sufficient  for  our  analysis  since  we 
are  observing  the  more  easily  manipulated  variables  in 


4b 


the  actual  design  of  a  particular  coating. 

Equation  119  is  readily  solved  numerically 
using  a  Newton  algorithm.  This  equation  is  solved  for 
the  N-layer  system  and  is  represented  on  a  three 
dimensional  plot.  The  axis  variables  are  layer  number 
I's  electrical  conductivity,  layer  number  I's  thickness 
and  differential  surface  temperature  increase.  Figures 
19  to  22  are  representative  plots  of  a  typical  2  layer 
system  at  various  incident  angles  and  polarizations. 

Later  in  the  experimental  verification  in 
Chapter  VI,  it  will  be  important  to  realize  that  a 
given  coating  will  be  considered  as  a  function  of 
incident  power  level  only.  The  coating  and  all  its 
particular  characteristics  are  already  fixed;  they 
cannot  be  varied  without  constructing  an  entirely  new 
coating  configuration.  Thus,  we  will  calculate 
equilibrium  temperature  as  a  function  of  incident 
power  level  using  equation  119  and  then  compare  these 
values  to  the  actual  measured  values.  Topologically 
then,  we  will  be  considering  single  lines  in  4N+7 


space . 
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COATINC  THICKNESS  IN  SKIN  DEPTHS 

PARALLEL  INCIDENCE:  Theta-  0  Degrees 


Numb..-  of  Lau.ri.  2  Pou.r  ir.  H-Cm  2-10  Fr,q.  ,n  CHi  =  2.45 
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Figure  19: 


Two  Layer  Electromagnetic  Heating  Profiles 
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1.0  3  0  4,0 


COflTJfW  THICKNESS  IN  SKIN  OtHTHS 

PARALLEL  INCIDENCE:  Theta“75  Degrees 


Nofftbff  of  tayers*  2  flower  in  W  Ciii'2  =  l0  Freq.  in  GHz“2.45 

Rir  Temp  •0*19.7  Convective  Exp*  .40  Emi  ss  i  •»!  tya  ,90  Plate  Height  'M>»  .100 
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Figure  20:  Two  Layer  Electromagnetic  Heating  Profiles 
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Hijint...'!  of  L»''e''f.=  i  Po'.i<i  It.  M  Cli.  ^=10  Fr^q.  in 
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Figure  21 i  Two  Layer  Electromagnetic  Heating  Profiles 
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Figure  22:  Two  Layer  Electromagnetic  Heating  Profiles 


51 


CHAPTER  VI 

EXPERIMENTAL  VERIFICATION 

The  coupling  of  the  electromagnetic  solution  with 
the  thermodynamic  solution  for  a  system  of  N  layers 
yields  a  theoretically  valuable  tool  for  predicting  the 
surface  temperatures  of  a  layered  system.  However, 
before  this  development  can  be  used  with  any  certainty, 
we  must  experimentally  verify  the  accuracy  of  the 
predicted  results.  The  overall  verification  was  a  three 
step  process.  We  first  measured  the  convective  heat 
transfer  exponent,  6,  and  then  examined  five  different 
conductive  samples  as  they  were  exposed  to  various 

2 

microwave  power  levels  ranging  from  10  to  40  mW/cm  . 
Irradiation  frequency  was  fixed  at  2.54  GHz.  An 
appropriate  convective  heat  transfer  coefficient  was  then 
chosen  that  would  allow  the  best  possible  fit  between 
experimental  and  theoretical  data. 

The  convective  heat  transfer  exponent  was 
measured  in  the  following  manner.  A  constant  power  was 
generated  within  a  thin  conductive  coating  (aquadaq)  by 
keeping  a  constant  direct  current  voltage  applied  accross 
the  sample.  This  sample  was  insulated  with  two  inches  of 
styrofoam  insulation  on  the  sides  and  back  so  that  we 


52 


could  closely  approximate  our  model  which  assumed  no 
conductive  heat  losses.  The  sample  was  allowed  to 
reach  steady  state  conditions  before  the  surface 
temperature  was  measured.  Surface  temperature  measure¬ 
ments  were  made  with  an  infrared  camera  capable  of 
measuring  temperature  variations  to  within  one  tenth 
degree  Centigrade.  We  then  calculated  an  absorbed  power 
density  given  by 

=  VI/A  (123) 

cL 

where 

V  =  applied  voltage  across  sample  (volts) 

I  =  current  through  the  sample  (amperes) 

A  =  area  of  the  sample  (meters  squared) 

This  absorbed  power  must  be  dissipated  via  surface 
convection  and  infrared  radiation;  therefore,  as  in  the 
development  in  Chapter  5,  we  equated  the  power  gains  to 
the  power  losses  which  resulted  in 


^  (“o-Ualr) 
H 


1  +  6  44 

"V(«o-«air) 


(124) 


Surface  emissivities  have  been  measured  by  other 
23  24 

researchers  ’  and  are  thus  known  quantities.  If  we 
assume  the  convective  coefficient  (1.42)  remains  relatively 
constant,  the  only  variable  remaining  is  the  convective 
exponent  which  may  be  calculated.  Its  value  is  given  by 
In 


“lirj) 


(125) 


Figure  23  illustrates  the  general  layout  of  the 
experiment • 
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Figure  23:  Schematic  of  Convective  Exponent  Experiment 


Test  equipment  used  in  the  experiment  consisted 
of  the  following: 

Nomenclature  Identification 

Power  Supply  Regulated  D.C.  Power  Supply 

0-50  VDC,  0-1.5A 
Kepco  Manufacturing  Co. 
Flushing,  New  York 

Digital  Multimeter 
Model  3466A 

Hewlett-Packard  Mfg.  Co. 
Colorado  Springs,  CO 


Ammeter 
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Voltmeter  Digital  Multimeter 

Model  3466A 

Hewlett-Packard  Mfg.  Co. 
Colorado  Springs,  CO 

Infrared  Camera  Thermovision  Model  680 

AGA  Manufacturing  Co. 
Secaucus,  New  Jersey 

The  sample  measured  consisted  of  a  10  cm  square 
of  conductive  coating  59  microns  in  thickness  and  mounted 
vertically.  Room  temperature  was  maintained  at  20.5f.l  C 
throughout  the  experiment.  Five  different  direct  current 
power  levels  were  observed  to  determine  the  linearity  of 
the  convection  exponent  with  temperature.  Figure  24 
illustrates  the  experimental  results  in  which  the  average 
value  of  the  convective  exponent  was  determined  to  be 
0.53.  Having  a  value  for  the  convection  exponent,  we 
then  verified  the  N-layer  electromagnetic  interaction. 
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There  are  two  major  shortfalls  of  the  N-]ayer 
model  which  have  a  drastic  impact  on  the  verification. 

The  first  of  these  is  the  semi -inf inite  plane  assumption. 
IVe  obviously  cannot  construct  a  semi -inf inite  layered 
system  in  the  laboratory  and  irradiate  it  uniformly 
with  a  plane  electromagnetic  wave;  therefore,  the 
approach  is  to  use  a  very  small  (much  less  than  a  wave¬ 
length)  multilayered  sample.  By  avoiding  resonant  shapes 
and  sizes  the  sample  will  experience  a  relatively  uniform 
electric  field  across  the  surface  and  the  observed 
surface  temperature  will  approximate  that  obse:^ved  on  a 
semi- inf inite  sheet. 

Secondly,  from  the  thermodynamics  point  of  view, 
we  assumed  the  conductive  coating  would  be  placed  on  a 
perfect  insulator  and  thus  no  thermal  conduction  would 
occur.  Obviously,  perfect  thermal  insulators  do  not 
exist  either;  hence,  we  have  another  source  of  error.  To 
minimize  the  thermal  conduction  losses,  samples  were 
placed  on  styrofoam  blocks  with  a  minimum  of  two  inches 
of  insulation  on  all  but  the  front  surface.  Styrofoam 
proved  to  be  well  suited  for  this  purpose  since  it  is  an 
excellent  insulator  and  has  a  measured  electrical 
permittivity  of  approximately  1.1  (10  GHz).  Electric¬ 
ally,  therefore,  the  foam  was  virtually  invisible  while 
thermally  it  provided  the  desired  insulation.  The 
samples  consisted  of  1,  2,  and  3  layer  configurations  all 


57 


cut  to  be  1.5  cm  square.  In  addition  to  the  equipment 
already  listed,  the  following  instruments  were  also 
used : 

Nomenclature  Identification 

Microwave  Generator  Microwave  Generator 

2.45  GHz,  0-200  Watt 

Kiva  Instrument  Co. 
Rockville,  Maryland 

Power  Meter  Radiation  Hazard  Meter 

(RAHAM),  Model  481 
General  Microwave 
Farmingdale,  New  York 

Figure  25  illustrates  the  physical  arrangement  of  the 
experiment . 


Insulator 


Figure  25:  Schematic  of  N-Layer  Verification 
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The  procedure  itself  was  relatively  straight 
forward.  The  microwave  power  meter  was  used  to  measure 
the  free  field  power  level  of  the  2.45  GHz  radiation. 
After  the  power  was  adjusted  to  the  proper  value  the 
layered  sample  was  placed  at  the  point  where  the  power 
was  measured.  The  samples  were  then  allowed  to  reach 
thermal  equilibrium  prior  to  measuring  their  surface 
temperatures  with  the  infrared  camera.  In  most  cases 
thermal  equilibrium  was  reached  in  less  than  10  minutes; 
however,  the  samples  were  irradiated  at  least  20  minutes 
each.  Each  sample  was  exposed  to  power  levels  of  10,  15, 
20,  30,  and  40  mW/cm  .  The  measured  temperature  for  each 
is  plotted  versus  input  power  in  figures  26  to  30. 
Calculated  values  are  also  plotted  on  the  same  graph  for 
comparison  purposes.  The  single  greatest  variation 
occurred  for  the  two  layer  aquadaq/plexiglas  combination. 
In  this  case  the  average  deviation  from  theoretical 
values  was  about  33  percent.  The  other  four  samples 
yielded  much  better  results  of  typically  less  than  20 
percent  variation.  To  obtain  the  theoretical  values 
shown  we  had  to  use  a  convective  coefficient  of  15  rather 
than  the  1.42  used  in  the  theoretical  discussion  in 
Chapter  V  thus  indicating  that  convection  is  much  more 
significant  than  originally  assumed. 


differential  surface  temperf 

( Cent i gr  ade ) 


Figure  26;  One  Layer  Electromagnetic  Heating 


60 


Figure  27:  One  Layer  Electromagnetic  Heating 


DIFFERENTIAL  SURFACE  TEMPERF 
(Cent (grade ) 


Figure  28;  One  Layer  Electromagnetic  Heating 


DIFFERENTIAL  SURFACE  TEMPERATURE 
(Cent (grade ) 
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MICROWAVE  INPUT  POWER  (Wsttz-Squara  Centimeter) 
(Frequency*'2 . 45  GHz) 


Figure  29:  Two  Layer  Electromagnetic  Heating 


DIFFERENTIAL  SURFACE  TEMPERATURE 
( Cent i grede ) 


MICROWAVE  INPUT  POWER  (Watt/Square  Centimeter) 
(Frequency*? . 45  GHz) 


Figure  30:  Three  Layer  Electromagnetic  Heating 
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CHAPTER  VII 

APPLICATION  OF  MDDEL  TO  TWO-DIMENSIONAL  TARGETS 

As  vrais  pointed  out  earlier  in  the  initial  development 
of  our  coating  model,  a  possible  limitation  to  its  use  results 
from  the  one -dimensional  analysis  itself.  This  chapter  discusses 
the  problem  of  applying  the  one- dimension!  coating  model  to  the 
more  "real  world"  two-dimensional  problem.  In  partici'lar,  the 
experimental  verification  itself  uses  an  electrically  small  two- 
dimensional  square  sanple  to  approximate  the  current  present  in 
the  semi-infinite  one- dimes ional  case.  Also,  the  question  of 
"nearest  neighbor"  influence  is  important  when  considering  elec¬ 
trically  large  sanples.  That  is,  v^en  we  observe  a  temperature 
at  a  point  cm  the  surface  of  a  coated  shape,  how  well  does  that 
relate  to  the  currents  present  directly  under  it?  Clearly,  the 
currents  in  the  neighboihood  of  this  point  must  contribute  some¬ 
thing  to  the  coating  temperature  but  the  question  is,  "how  signif¬ 
icant  is  their  influatice?"  In  order  to  qualitatively  address 
each  of  these  situations,  we  will  first  consider  the  electri¬ 
cally  small  sample  as  applied  to  the  verification  experiment 
and  later  discniss  the  "nearest  neighbor"  considerations  in 
regards  to  electrically  large  models. 

As  was  indicated  earlier  in  Chapter  VI,  the  verification 
process  for  the  semi -infinite  one -dimensional  model  is  a  difficult 
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one.  Physically,  we  were  constrained  in  size  by  an  anechoic 
chamber  with  a  maximum  usable  width  of  approximately  one  meter. 

In  the  far  field  we  had  an  area  of  approximately  20  x  15  cm  in 
which  the  microwave  power  density  was  relatively  uniform.  Elec- 
trically  we  could  achieve  a  power  density  of  20  "W/cm  in  the 
above  area  at  a  frequency  of  2.45  (Hz.  The  radiation  pattern  was 
formed  by  use  of  a  quarter  wave  monopole  in  conjunction  with  a 
parabolic  reflector  over  an  aluninum  ground  plane.  See  figures 
31,  32,  and  33  below. 


Figure  31:  Photograph  of  Monopole  and  Parabolic  Reflector 
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Figure  33:  Anechoic  Chanter  With  Electrically  Small  Target  and 
Infrared  Camera  in  Foreground 
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The  verification  problem  is  such  that  it  could  proceed 
in  either  of  two  ways.  The  first  was  to  construct  several  elec¬ 
trically  very  large  sanples  that  by  nature  of  their  size  would 
the  smi -infinite  condition  of  our  model.  There  are 
several  problems  associated  with,  this  approach.  Large  coated 
targets  are  considerably  difficult  to  construct  since  they  must 
have  completely  uniform  conductive  coatings  over  their  surfaces. 

The  problem  is  even  more  difficult  with  two  or  three  composite 
layers.  The  largest  problem,  however,  results  from  the  nature  of 
the  anechoic  chamber  itself.  That  is,  if  we  placed  a  large 
target  completely  across  the  diamber,  nearly  all  of  the  incident 
microwave  field  would  be  back  scattered  into  the  parabolic 
onitter  vMch  in  effect  would  turn  our  anechoic  chamber  into 
some  type  of  tuned  cavity.  Thus,  we  would  know  nothing  about  the 
fields  present  at  the  coating.  The  back  scatter  problem  could  be 
reduced  or  eliminated  by  removing  the  parabolic  reflector  and 
substituting  anechoic  absorbers  but  then  our  incident  field 
strength  would  have  suffer^  drastically.  Therefore,  because  of 
the  mechanical  and  electrical  difficulties  inherent  with  the  large 
model  approximation,  it  appears  more  desirable  to  use  an  elec¬ 
trically  small  model  to  approximate  the  send -infinite  results. 

Ihe  electrically  small  sample  has  several  advantages. 
First,  it  is  much  easier  to  fabricate  and  characterize  elec¬ 
trically.  In  fact  we  can  measure  the  coating  thickness  directly 
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using  a  good  micrometer  since  the  micrometer  jaws  nearly  com¬ 
pletely  cover  the  sample's  face.  The  back  scatter  problem  is 
practically  nonexistent.  Similar  to  developments  elsewhere,  if 
we  consider  a  sample  vdiich  is  electrically  small  and  assume  the 
incident  field  induces  a  relatively  uniform  surface  current 

across  the  sample  faces  then  we  may  treat  it  as  a  Huygen's  source 

25 

(We  will  discuss  the  uniformity  assumption  more  later.)  .  See 
figure  34  below. 


Figure  34:  A  IMiform  Current  Sheet  Radiator 
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If  we  further  assume  the  induced  current  is  not  only  uniform  but 
also  the  same  as  that  induced  on  an  infinite  conductor  and  that 


contour  charge  is  negligible,  then  we  have  The  radia- 


tion  field 

is  *  -joiX 

(126) 

where 

pJ„dy)  d  X  . 

(127) 

b:*- 

jo)  cosiji  cose 

(128) 

Eg  =  j 

u)  Ajj  sin<|) 

(129) 

Therefore,  the  approximate  magnitude  of  the  backscattered  field 
on  the  z  axis  and  polarized  in  the  x  direction  (9*0,  (}>  *  0) 
would  be 

E .  dy  dx 
l^slo  “  “  ^-Jyrc-T-^ 


Asstaning  that  this  field  is  redirected  by  the  parabolic  reflector 
with  some  gain  G,  the  redirected  field  strength  at  the  target 
would  be 


|Etl  ■ 


u)GE.  dy  dx 
8irc  r 


(131) 


The  percentage  of  re-reflected  to  incident  power  at  the  target 
face  would  be 


(132) 


For 
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w  =  2it  X  2.45  X  10^  Hz 

dy  =  dx  =  .015  m 
G  =  20 
r  =  1.3  m 
c  =  3  X  10®  m/s 
[E  [  ^ 

we  have  — a,  5  x  lO'^  %  (133) 

|Eil 

In  other  words,  our  small  sample  size  interferes  with  the  incident 
field  only  a  small  amount.  One  obvious  advantage  of 
this  is  that  we  can  measure  and  characterize  the  microwave  fields 
within  the  anechoic  chamber  with  no  targets  present  and  then  be 
sure  that  this  characterization  changes  very  little  with  the 
addition  of  a  small  object. 

Two  other  considerations  that  must  be  addressed  before 

the  small  sample  may  be  used  with  confidence  concern  the  surface 

current  uniformity  and  magnitude  on  such  targets.  Considerable 

insight  is  provided  to  both  of  these  questions  in  an  article  by 

Wilton,  et  al,  in  which  a  method  of  moment  solution  is  considered 

for  various  size  flat  conducting  plates.  See  figures  35  and 

36  below  for  a  profile  of  surface  currents  on  a  l.OX  and 
26 

0.15X  square  plate'^  . 
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The  ijjpjrtant  aspects  to  note  are  the  magnitude  of  J  as  the  plate 

size  decreases  and  also  the  more  uniform  current  distribution 

as  the  plate  size  reduces.  In  an  e3q>eriment  to  observe  these 

effects  more  graphically,  l.OA,  O.Sx,  and  0.25x  square  plates 

were  irradiated  with  normally  incident  microwaves  at  a  power 
2 

level  of  20  mV^cm  .  See  figures  37  to  45. 


Figure  37:  Photograph  of  l.OX  Square  Plate 
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Figure  38:  Infrared  Photograph  of  a  l.OX  Square  Plate 


Figure  39:  Thermal  Profile  Across  the  G?nter  of  a  l.Ox 


Square  Plate 


Figure  40:  Actual  Photograph  of  a  O.SX  Square  Plate 


Figure  41:  Infrared  Photograph  of  a  0.5X  Square  Plate 


Figure  42:  Thermal  Profile  Across  the  Center  of  a  O.Sx 


Square  Plate 


Figure  43:  Actual  Photograph  of  a  0.25X  Square  Plate 


Figure  44:  Infrared  Photograph  of  a  0.25\  Square  Plate 


Figure  45:  IheTinal  Pixfile  Across  the  Center  of  a  0.25X 
Square  Plate 
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The  ligjiter  areas  correspond  to  greater  surface  currents. 
A  single  temperature  profile  was  taken  across  the  middle  of  each 
target  to  allow  comparison  with  the  method  of  moments  results  by 
Wilton.  Notice  the  great  similarity  between  the  l.OA  e:q)eri- 
mental  plate  and  theoretical  results.  To  consider  the  question 
of  uniformity  more  graphically.  Figure  46  is  a  plot  of  a  .ISx 
square  plate  in  which  the  shading  depicts  those  areas  in  which 
the  surface  current  varied  less  than  20%  from  that  vrtiich  would 
be  observed  on  an  infinite  conducting  sheet.  This  information 
is  graphically  taken  from  the  plot  of  Wilton,  Figure  36. 


I  p 

•  incident 
Shading:  ±  20% 


Figure  46:  Current  Uniformity  on  a  .15x  Plate 
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Therefore,  by  observing  the  ten^erature  at  the  center  of  our 

small  sample,  we  are  reasonably  assured  that  our  results  are 

similar  to  vAraX  would  be  observed  on  a  semi -infinite  plate. 

In  the  experimental  verification  snail  square  saii;>les  were 

used  which  facilitated  the  above  discussion;  however,  this  is  not 

27 

an  absolute  requirement.  Andre jewski  has  confuted  the  surface 

current  at  the  center  of  a  disc  as  a  function  of  c  where 
c  =  with  "a"  being  the  radius  of  a  sample  disc.  See  figure 

47  for  a  plot  of  his  results. 


Figure  47:  Current  at  the  Center  of  a  Disc  as  a  Function  of  c 
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As  in  the  case  of  the  sn\all  square,  notice  that  as  c  gets  small 

the  field  at  the  center  of  the  disc  approaches  the  semi-infinite 

plane  value  (  li?|=l  ).  In  fact  for  a  disc  1.5  cm  in  diameter  the 

predicted  surface  current  value  at  the  center  is  within  10%  of  the 

value  experienced  by  a  semi-infinite  conducting  sheet.  Therefore, 

it  seems  reasonable  that  we  may  use  small  discs  as  well  as  squares 

in  the  verification  process  with  no  loss  in  validity.  For  a  more 

con^lete  discussion  of  theoretical  solutions  of  the  disc  one  may 
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refer  to  any  of  several  good  references.  *  * 

All  the  discussion  to  this  point  has  been  limited  to 

the  small  sanples  used  in  the  verification  process.  We  have 
developed  a  model  and  have  shown  that  it  is  valid  in  regards  to 
whether  or  not  a  particular  coating  configuration  will  heat 
properly  when  exposed  to  microwave  radiation.  One  final  question 
vfliich  is  relevant  vdien  considering  electrically  large  samples 
is  the  question  of  "nearest  neighbor"  influence.  That  is,  how 
well  can  our  coating  model  predict  surface  currents  below  it? 

Recall  that  in  general  our  coating  will  be  displaced 
from  the  surface  of  the  object  we  are  observing.  Although,  this 
displacement  is  necessary  in  order  to  observe  heating,  it  is 
an  overall  detriment  to  the  resolving  power  of  the  coating. 

We  may  begin  to  understand  th;  s  problem  more  if  we  recognize 
that  in  order  for  a  parti  oilai  current  element  to  interact  with 
the  coating,  that  current  elaient  mist  emit  an  electromagnetic 
wave  which  in  turn  must  trave.  to  the  coating.  See  Figure  48. 


Neighbor  Contribution 

Desired  Current  Measurement 
.  SURFACE  CURRENT 


*  (Out  of  Page) 


INSULATING  LAYER 


Measurement  Point 
(Infrared  Radiation) 


Figure  48:  Cross  Section  of  a  Conductive  Coating  Displaced  Away 
From  Condhictor  Surface 

Ihe  question  then  is  hovr  far  away  from  the  measurement  point 

tan  this  neighbor  contribution  be?  We  may  receive  some  insight 
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ix>  the  problffln  by  considering  the  typical  waveguide  solution 
to  two  semi -infinite  planes  separated  by  a  distance  "a”.  See 
figure  49  below. 
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Assuning  wave  propagation  in  the 

solution  is  given  as 

z  direction  and  a  TE  mode,  the 

Ey  =  sin  x)  e'"^^ 

(134) 

m  =  1,2,3,  . . . 

(135) 

(136) 

We  see  that  y  vdll  either  be  pure  imaginary  or  real  depending 
on  the  argmient  under  the  radical.  If  y  is  real,  the  wave  will 
attenuate  exponentially  vdiich  is  vdiat  we  desire  in  order  to 
reduce  the  nearest  neighbor  effect.  We  may  solve  for  a  minimum 
value  of  a  in  which  this  will  occur.  Forcing  the  argument  to  be 
greater  than  zero  and  solving  for  a,  we  find  that  attenuation 
will  occur  as  long  as 

a  <  — S —  (137) 

If  we  assune  the  sin^lest  propagation  mode  Cm  ^  1)  and  a  foam 

insulating  layer  (e^  =  Wj.  “  1) »  then  we  have  lateral  attenuation 

as  long  as  the  insulator  is  less  than  thick.  Obviously  the 

attenuation  will  be  stronger  as  the  layer  gets  thinner.  The 

distance  at  which  an  attenuated  wave  has  reductsd  to  a  value  ^ 

of  its  initial  value  is  given  by 
1 _ 

-  (D  we 


z 


(138) 
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In  the  regime  where  we  nwst  operate  in  order  to  have  significant 
lateral  attenuation  and  thus  little  nearest  neighbor  influence, 
we  have 


(m.,2  2 


»  0)  pe 


(139) 


Therefore, 

z  =  ^  (140) 

mir 

In  other  words,  we  would  not  expect  to  see  significant  nearest 
neighbor  influence  at  distances  much  over  a/ir  for  the  worse  case 
of  m=l. 


An  experiment  was  conducted  in  order  to  further  understand 
the  significance  of  the  nearest  neighbor  problems.  A  large 
aluminum  conductor  was  covered  with  a  layer  of  foam  5.1  ram  in 
thickness.  The  foam  was  subsequently  coated  with  an  aquadaq 
coating  approximately  15  microns  in  thickness  and  electrical 
conductivity  of  approximately  315  mhos/m.  As  the  large  target 
was  irradiated,  1.5  cm  diameter  holes  were  placed  in  the  aluminum 
at  closer  and  closer  intervals.  The  holes  behind  the  coating 
were  observed  as  hot  spots.  A  point  was  finally  reached  at  which 
time  the  two  spots  had  begun  to  merge.  At  this  point  the  centers 
of  the  holes  were  2.0  cm  apart.  See  figures  50  to  61  for 
thermovision  photographs  and  thennal  profiles  of  the  results. 
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Figure  50:  Photograph  of  Electrically  Large  Plate  in  Place  in 
the  Anechoic  Qiamber 


Figure  51:  Photograph  of  the  Single  1.5  an  Fble  in  the 
Alimimin  Plate 


Figure  52:  Photograph  of  the  Dcnible  Hole  Configuration 


Figure  53:  Infrared  Photograph  of  Artificially  Heated  Double  Hole 
Configuration  for  Size  Cotnparison  with  Microwave  Results 


Figure  54:  iTifrared  Photograph  of  Single  Hole  Heating  Pattern 
Resulting  from  Microwave  Radiation  (1°  C  Scale) 


Figure  55:  Horizontal  Thermal  Profile  of  Single  Hole  Heating 
Pattern  (1°  C  Scale) 
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Figure  58:  Infrared  Results  of  Double  Hole  Microwave  Heating 
Pattern  (1°  C  Scale) 


Figure  59:  Horizontal  Thermal  Profile  of  Double  Hole  Heating 
(1°  C  Scale) 
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Figure  60; 


Same  as  Figure  58  Except  2°  C  Scale 


Figure  61: 


Same  as  Figure  59  Except  2°  C  Scale 
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A  z  calculated  for  the  siii^slest  On®l)  mode  in  this 
configuration  is  approximately  .2  cm.  Note  that  we  may  still 
observe  separation  oi  the  two  holes  on  the  themovision  photo¬ 
graphs  and  at  that  point  the  actual  metal  separating  the  holes 
was  .5  cm.  It  is  of  note  also,  that  the  particular  infrared 
camera/object  separation  distance  for  this  e3q)eriment  only  provided 
a  .17  cm  camera  resolution. 

As  a  final  note,  it  should  be  obvious  that  electrically 
large,  finite  dimensional,  targets  exhibit  large  variations  in 
surface  current  density  (See  figures  35,38,  and  39  for  examples,). 
It  is  possible  to  gain  some  insight  into  how  these  variations 
relate  to  the  one-dimensional  model  by  considering  the  following. 
Suppose  that  we  have  a  two-dimensional  problem  in  which  there  are 
surface  current  variations  in  the  x  direction  of  a  substrate  as 
a  result  of  its  finite  size.  See  Figure  62. 
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Qearly,  if  the  substrate  were  infinite  in  the  x-y  plane, 
the  value  of  J  would  be  a  constant  given  by  .  However,  as  a 
result  of  its  finite  size  in  the  x  direction  and  the  inqjosition 
of  boimdary  conditions,  we  actually  have  a  current  distribution 
varying  spatially  in  the  x  direction  and  si;5)posed  to  be  that 
illustrated  in  Figure  62.  In  this  example  only,  the  surface 
current  is  confined  to  the  range  0<J<3Hj^. 

We  know  that  the  total  electric  field  in  the  coating 
will  be  the  vector  sum  of  the  incident  field  plus  the  scattered 
field  from  the  current  element  beneath  the  measurement  point,  P. 
The  question  of  nearest  neighbor  contribution  to  the  field  at  P 
has  already  been  discussed.  It  was  shown  that  most  of  the 
scattered  field  at  P  resulted  from  current  distributions  within  a 
circular  radius  of  approximately  a/ir.  Therefore,  if  a  is  small, 
we  may  assume  the  radiating  source  to  be  located  directly  behind 
P.  We  have 

=  Si  *  63  (141) 

vdiere  is  the  field  in  the  coating,  is  the  incident  field, 

and  Sg  is  the  scattered  field.  We  may  define  a  fourth  field,  2^, 

as  the  scattered  field  from  an  infinite  substrate  and  2^  as  a 

P 

perturbing  field  that  when  combined  with  2^^  would  yield  the 
actual  scattered  field  from  the  finite  size  substrate.  Thus, 

S3  -  S,  .  ^  (142) 

Substituting  equation  142  into  equation  141  we  have, 

C  1  “  p 

Recalling  the  development  in  Chapter  V,  we  may  write  an  expression 
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for  the  deposited  power  density  in  the  coating,  as  follows: 

'5c  ■  ■  4  ■ 

Taking  the  dot  product  of  with  itself,  we  have 

«  a{  ^  +  2(  )}  (145) 

Rearranging,  we  have 

»  o{  eJ  +  +  2\\  ^  C  E^  +  2Ct-tp<-^p))}  (146) 

Qearly,  by  considering  tenperatijre  regimes  in  which  no 
changes  of  state  occur  for  our  coating  and  by  remembering  the 
conservation  of  energy,  we  may  infer  that  the  tenqjerature  of  the 
coating  is  a  continuous,  increasing  function  of  Q^;  that  is, 

4T  =  f(Q^.,H)  (147) 

where  H  consists  of  several  thermodynamic  variables  such  as 
surface  emissivity,  san^le  height,  orimtation,  roughness,  etc.. 

The  ijnwrtant  point,  however,  is  that  for  each  increment  we 
Increase  Q^.,  AT  will  also  increase.  If  we  consider  a  simple 
case  in  vfluch  the  power  loss,  Qj^,  may  be  expressed  in  terms  of 
Newton's  law  of  cooling,  we  have 

Ql  -  gCH)AT  (148) 

’jdiere  g(H)  is  some  function  of  the  thermodynamic  variables 
listed  above;  for  a  given  coating/substrate  configuration 
it  is  assimed  that  g(H)  will  be  relatively  constant  independent 
of  tei^>erature.  For  exan^jle,  infrared  emissivity  is  a  function 
of  coating  electrical  corxluc tivity;  electrical  conductivity 
was  measured  vers^..  o  ^so'  ..e  temperature  for  an  aquadaq  coating 
in  Appendix  A  and  found  to  vary  less-  than  0.3%  over  the  temperature 
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ranges  of  interest.  It  can  be  shown  that  the  anissivity  would 
vary  even  less  than  the  conductivity;  thus,  it  is  assumed  to  be 
constant.  We  may  now  write  the  conservation  of  energy  equation 
for  the  coating/substrate  system  in  dynamic  equilibrium;  we  have 
Ql  -  (149) 

vdiere  t  is  the  coating  thickness.  In  equation  149  we  have 
assumed  that  the  coating  in  this  example  is  so  thin  that  the  power 
deposition  is  uniform  throughout  and  is  only  a  function  of  the 
coating  separation  from  the  substrate.  Substituting  equation 
148  into  equation  149,  we  find 

4T  “  Qj.(z“-a)  (150) 

This  is  the  sou^t  after  result;  the  differential  surface 
tenperature  is  directly  proportional  to  the  deposited  power 
density.  The  importance  of  this  arises  from  the  composition 
of  (^.  The  first  part  of  equation  146,  (  E?  )» 

is  vdiat  we  might  call  the  first  order  terms  vMch  result  from 
a  purely  one-dimensional  analysis  of  the  problem  (  The  coating/ 
substrate  cooibination  is  infinite  in  the  x-y  plane.).  The 
remainder  of  equation  146,  (  ^  +  ^^’(  ^ 

interpreted  as  the  higher  order  contributions  resulting  from 
the  imposed  boundary  conditions  in  a  finite  case.  Thus,  the 
differential  ten^jeratures  we  observe  on  our  coating  surface 
may  be  thought  of  as  the  sum  of  tenqieratures  we  would  observe 
on  an  infinite  configuration  plus  any  changes  that  would  occur 
as  a  result  of  observing  a  finite  size;  that  is, 


at  ■  AT  +  aT_ 

CO  p 


(151) 
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In  equation  151  bear  in  mind  that  aT  may  have  a  negative  value 
since  only  AT  is  required  to  be  greater  than  zero  by  the  law 
of  conservation  of  energy.  AT^  is  the  term  that  results  from  the 
interference  of  our  pertinrbing  field,  with  the  total  field 
that  would  be  present  in  the  case  of  an  infinite  coating/ substrate 
configuration,  C  Also,  AT  typically  is  not  directly 

proportional  to  the  absorbed  power  density  but  rather  approximates 
a  square  law  to  a  large  degree.  This  non-linearity  is  no  problem 
since  the  infrared  data  is  generally  entered  directly  into  a 
computer  and  may,  therefore,  be  corrected  for  in  a  relatively 
straight  forward  manner. 

The  iii5X)rtance  of  all  this  is  that  in  analyzing  the 
current  distribution  on  a  finite  plate  via  the  coating 
ten^jerature  distribution  we  may  wish  to  consider  only  the  non- 
steady  state  term,  AT^.  Ihat  is,  in  the  process  of  completing  an 
infrared  measurement  we  may  wish  to  scale  the  results  by  some 
additive  constant  in  order  to  correlate  the  infrared  results 
with  a  single  probe  measurement  on  a  particular  surface  and 
thereby  "calibrate"  the  infrared  results  for  all  other  points 
on  that  surface. 

Since  we  have  shown  earlier  that  small  samples  may  be 
used  to  approximate  the  one  dimensional  infinite  plane,  it 
seems  reasonable  then  to  assume  that  we  can  use  small  samples 
to  measvire  AT^  directly.  A  series  of  flat  plate  experiments  were 
performed  in  which  this  was  done.  First,  a  small  1.5  cm  square 
was  irradiated  with  microwaves  at  a  frequency  of  2.45  GHz  and  a 
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power  of  25  mWcm  .  This  is  the  sanse  sized  san5>le  iised  in  the 
verification  experiments,  TTierefore,  by  measuring  the  equilibrium 
temperatures  on  the  surface  of  the  mnall  sanple  we  have  obtained 
for  tfiat  particular  coating  configuration.  Next,  a  l.Ox 
square  plate  was  irradiated  in  order  to  illustrate  the  process. 

An  infrared  photograph  was  made  along  with  a  single  thermal 
profile.  ATp  may  be  obtained  by  subtracting  from  the  total 
temperature  measured  by  the  thermovision.  AT^  was  drawn  on  the 
profile  photograph  to  illustrate  the  magnitude  involved.  Figures 
63  to  68  were  made  from  a  coating  of  15  microns  of  aquadaq 
(315  mhos/m)  placed  on  5.1  mm  of  foam  insulator.  Figures  69 
through  74  were  made  using  an  aquadaq  coating  15  microns  thick 
and  placed  on  10.2  mm  of  foam  insulation.  Additionally,  Figures 
69  through  74  included  a  copper  substrate  whereas  the  earlier 
photographs  contained  no  substrate. 


Figure  63:  Thermovision  Photograph  of  1.5  cm  Square  Sanqile 
CNo  Substrate) 


Figure  64:  Thermal  Profile  (Horizontal)  through  the  Center  of 
the  1,5  cm  Square  San5)le  (2,45  GHz  at  25  mw/cm  ) 


Figure  65:  Ihernovision  Photograph  of  l.OV  Square  Plate  (Horizontal 
nMte  line  indicates  area  Where  thennal  profile  was 
taken.) 


Figure  66:  ITiennal  Profile  from  l.O^  Square  Plate  with  AT, 
Plotted 
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Figure  67;  Thermovision  Photograph  of  O.SX  Square  Plate  (Horizontal 
white  line  indicates  area  where  thermal  profile  was 
taken.) 


Figure  68;  Thermal  Profile  from  0,5X  Square  Plate  with  AT 
Plotted 


Figure  69:  Themovlsion  Ffiotognqph  of  1.5  on  Square  Sample 
(Copper  Substrate) 


Figure  70:  ‘Hieniial  Profile  G^wizontal)  thorou^  the  Center  of 
the  1.5  cm  Scpiare  Sanple  (2.45  GHz  at  25  mt/car) 
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Figure  71:  Ihenaovisioii  Riotograph  of  l.OX  Square  Plate  with 
a  Copper  Substrate 


Figure  72:  Ihennal  Profile  Corresponding  to  White  line  in 
Figure  71 


Figure  73;  Thermovision  Photograph  of  a  0.5\  Square  Plate  with 
a  Copper  Substrate 


Figure  74:  Thermal  Profile  Corresponding  to  White  Line  in 
Figure  73 
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In  passing  it  should  be  pointed  out  that  the  preceding 

results  are  consistent  with  the  theoretical  analysis  typically 

employed  in  the  development  of  holography  ,  This  is  to  be 

expected  since  the  coating  is  siu5)ly  recording  the  interference 

patterns  resulting  from  the  combined  effects  of  the  incident  and 

scattered  fields  from  the  substrate  surface.  This  would  be  the 

37 

arrangement  for  a  typical  Gabor  hologram  .  As  is  the  case  of  an 

optical  hologram,  we  have  both  phase  and  anqplitude  information 

stored  in  the  thermogram  of  the  surface  coating.  Qearly,  the 

points  on  a  thermal  profile  lying  below  AT^  indicate  that  tp 

must  be  out  of  phase  with  which  would  yield  a  negative 

value  for  the  dot  product  in  equation  146.  It  is  not  as  easy 

to  assess  phase  information  for  the  aT’s  above  aT^,  These  areas 

may  correspond  to  ^p’s  out  of  phase  with  but  much  greater 

in  magnitude  such  that  ?p  itself  is  the  primary  heating  source. 

Continued  research  in  this  area  may  finally  provide  much  needed 

information  in  regards  to  surface  current  phase.  Fast  Fourier 

Transform  (FFT]  techniques  and  optical  reduction  of  the  coating 

interference  patterns  into  optical  wavelengths  for  subsequent 

laser  illumination  are  both  techniques  that  may  provide  the  desired 
38 

phase  information. 

This  was  not  an  attempt  to  quantitatively  determine 
the  resolving  power  of  any  particular  coating.  That  effort 
is  better  left  as  a  topic  for  further  study  on  its  own  merit 
in  another  dissertation.  However,  it  does  illustrate  that  there 
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are  no  far  reaching  ’’nearest  neighbor”  effects.  We  observed  no 
shadows  or  false  apertures  as  might  be  expected  if  there  were 
some  other  propagation  mode  present.  Thus,  the  end  result  of 
this  series  of  experiments  is  a  qualitative  understanding  of  the 
nearest  neighbor  interaction  in  electrically  large  targets  plus 
an  assurance  that  the  one  dimensional  model  may  be  used  with 
confidence  in  the  design  of  a  particular  coating. 
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CHAPTER  VIII 

RECOMMENDATIONS  AND  CONCLUSIONS 

This  paper  provides  a  method  of  selecting  an 
optimum  coating  for  use  in  the  thermographic  detection  of 
microwave  induced  surface  currents.  Optimum  in  this  sense 
implies  two  things.  First,  we  needed  a  coating  that  would 
heat  to  at  least  1  degree  Centigrade  above  ambient 
temperature  to  achieve  10  levels  of  resolution  on  the 
infrared  system.  This  is  directed  by  the  accuracy  of  the 
thermographic  camera  which  can  resolve  temperatures  to 
within  0.1°  C.  Secondly,  we  needed  a  coating  that  would 
heat  to  an  acceptable  level  but  minimally  interfere  with 
the  fields  present.  This  second  requirement  stems  from 
the  necessity  of  minimizing  the  effect  any  measurement 
device  would  have  with  the  quantity  being  measured. 

In  attaining  an  optimum  coating  there  were  other 
constraints  that  we  had  to  adhere  to  in  order  to  have  a 
realizable  system.  First,  the  coating  should  not  contain 
ferrites,  non- isotropic ,  or  permeable  materials.  This 
requirement  was  established  to  ensure  that  the  coating 
itself  did  not  introduce  possible  intermodulation 
distortions  resulting  from  non-linear  materials. 


107 


Secondly,  it  was  desired  that  the  coating  be  easily 
fabricated  for  application  to  a  variety  of  possible 
shapes.  And  thirdly,  it  had  to  be  non-toxic  to  avoid 
possible  health  problems. 

From  an  engineering  point  of  view,  there  are  only 
a  finite  number  of  variables  that  may  be  controlled  to 
achieve  the  desired  results.  The  material  character^ 
istics  are  limited  to  coating  electrical  conductivity  and 
the  insulating  layer  permittivity.  Since  we  do  not  allow 
the  use  of  ferrous  materials  in  either,  the  relative 
permeability  is  approximately  1.  The  easiest  factor  to 
control,  and  the  most  effective,  is  the  individual  layer 
thickness . 

Before  explaining  the  coating  development,  it  is 
appropriate  to  point  out  the  two  different  uses  of  the 
coatings.  The  first  involves  a  situation  in  which  only 
the  shape  of  a  particular  object  is  important  in  an 
electromagnetic  interaction  problem.  To  determine  the 
current  distribution  in  this  problem,  the  simplest  and 
most  effective  technique  is  to  construct  a  model  of  the 
object  from  a  foam  material  (low  thermal  mass)  and  then 
coat  it  with  a  thick  (2-3  skin  depths)  coating  of  high 
electrical  conductivity  (300-500  mhos/m)  material.  Thus, 
we  simulate  a  highly  conductive  object  which  may  be  made 
of  something  such  as  aluminum  with  one  of  lesser 

r  2 

conductivity.  Senior,  et  al 


have  demonstrated 
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that  the  modeling  would  not  be  one-to-one  for  small 
conductivities  but  Sega  has  demonstrated  that  for  reason¬ 
able  conductivities  (greater  that  300  mhos/m)  the  corapari- 

V  ,  33  . 

son  may  be  close  .  A  possible  application  of  such  a 
scheme  might  be  to  determine  the  current  distributions  on 
an  aircraft,  ship,  or  spacecraft  to  facilitate  the  proper 
placement  of  antennas  or  measure  radar  cross  section. 
Additionally,  this  type  information  might  be  valuable 
when  studying  the  vulnerability  of  an  existing  or  future 
weapon  system  to  the  effects  of  an  Electromagnetic  Pulse. 

The  more  difficult  problem  is  one  which  cannot  be 
modeled  because  the  materials  the  object  is  made  from  are 
relevant  factors  in  the  particular  microwave  interaction. 

An  example  is  a  missile  or  satellite  constructed  from 
conductive  composites,  iron  alloys,  aluminum,  etc.  Here 
the  object  itself  must  be  irradiated  and  the  resulting 
surface  currents  measured.  To  measure  these  currents,  a 
different  type  coating  arrangement  is  required.  Two 
requirements  exist;  the  coating  must  be  thermally  isolated 
from  the  surface  to  allow  heating,  plus  it  should  be 
sufficiently  distant  from  the  surface  so  that  the  total 
(reflected  plus  incident)  electric  field  is  reasonably 
large.  That  is,  the  electric  field  must  be  large  enough 
so  that  the  heating  in  the  coating  will  result  in  a 
temperature  increase  greater  than  one  degree  Centigrade. 
Neoprene  rubber  fills  several  of  the  requirements  for  the 
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insulating  layer.  It  is  a  good  thermal  insulator;  it  is 


non-ferrous;  it  will  easily  conform  to  complex  shapes 
with  the  aid  of  contact  cement;  and  it  has  a  high 
relative  permittivity  (measured  to  be  31  at  10  GHz) . 

The  high  permittivity  allows  us  to  place  the  coating  an 
electrically  long  distance  from  the  surface  even  though 
the  rubber  may  be  physically  thin.  The  advantage  of  this 
may  be  seen  if  one  considers  the  infinite,  perfectly 
conducting  plane.  At  the  surface  we  have  a  surface 
current  magnitude  given  by  2Hq,  but  the  electric  field  is 
zero  there;  however,  at  a  distance  of  one  quarter  wave¬ 
length  from  the  plane  the  electric  field  magnitude  is 
given  by  2Eq.  Thus,  if  we  placed  a  thin  conductive  layer 
at  this  distance  from  our  infinite  conducting  plane,  we 
would  see  maximum  heating.  A  large  permittivity  material 
allows  us  to  minimize  the  actual  insulator  thickness  thus 
making  it  easier  to  apply.  There  is  a  disadvantage,  how¬ 
ever,  in  that  the  high  permittivity  insulator  has  a  much 
higher  reflection  coefficient  which  has  the  effect  of 
shielding  the  surface  from  the  incident  microwaves. 

To  optimize  the  above  situation,  that  of  placing 
our  conductive  coating  on  an  insulating  layer,  we  may  use 
the  following  procedure.  A  computer  program,  Uthick, 
described  in  Appendix  G  ,  generates  a  square  matrix  in 
which  the  coating  thickness  and  insulator  thickness  have 
been  allowed  to  vary  for  the  corresponding  matrix 
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elements.  The  coating  thickness  changes  linearly  along 
the  rows  and  the  insulator  thickness  changes  linearly 
along  the  columns.  After  this  matrix  is  generated  for  a 
particular  conductivity  and  insulator  configuration,  the 
results  are  plotted  on  a  contour  plot.  The  contour  pro¬ 
gram,  Contor,  is  located  in  the  appendix  also.  Each  coat¬ 
ing  configuration  is  examined  with  and  without  a  conductor 
behind  the  insulating  layer,  thus  simulating  the  situation 
of  having  currents  present  and  not  present.  Finally,  for 
comparison  purposes  a  different  contour  is  presented  which 
indicates  the  thermrl  differences  between  a  coating 
scheme  with  and  without  the  conductor  behind  the 
insulator.  Therefore,  we  can  determine  at  a  glance  what 
the  optimum  coating  configuration  would  be  for  a  given 
situation.  For  electrical  conductivities  ranging  from 
0.1  to  40  rahos/m  the  appendix  on  electrical  conductivity 
describes  in  detail  how  a  particular  conductivity  coating 
may  be  made  using  a  carbon  and  paraffin  mixture.  Figures 
75  to  95  illustrate  the  various  contour  plots  for 
conductivities  of  1,  5,  and  10  mhos/m  with  insulators  of 
styrofoam  (permittivity  =  1.1)  and  neoprene  (permittivity 
=  31) .  Additionally,  a  plot  of  aquadaq  (conductivity  = 

315  mhos/m)  on  1/16  inch  plexiglas  is  given.  The  incident 

2 

microwave  field  was  set  at  10  raW/cm  and  2.45  GHz. 
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Figure  75 ;  Carbon/Paraffin/Styrofoam  Thermal  Contours 
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Figure  76:  Carbon/Paraf f in/Styrofoam/Copper  Thermal 
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Figure  77:  Carbon/Paraf fin/Styrofoam/Copper  Difference 
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Figure  78:  Carbon/Paraffin/Styrofoam/Thermal  Contours 
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Figure  85;  Carbon/Paraffin/Neoprene/Copper  Thermal 
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Figure  86:  Carbon/Paraffin/Neoprene/Copper  Difference 
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Figure  87 :  Carbon/Paraffin/Neoprene  Thermal  Contours 
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Figure  88  :  Carbon/Paraffin/Neoprene/Copper  Thermal 
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Figure  89:  Carbon/Paraffin/Neoprene/Copper  Difference 
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Figure  90:  Carbon/Paraffin/Neoprene/Thermal 
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Figure  92:  Carbon/Paraffin/Neoprene/Copper  Difference 
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Figure  93:  Aquadaq/Plexlglas  Thermal  Contours 
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Figure  94:  Aquadaq/Plexiglas/Copper  Thermal  Contours 
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Figure  95 


Aquadaq/Plexiglas/Copper  Difference  Contours 
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The  preceding  figures  were  included  to  primarily 
illustrate  the  effects  of  different  coating/insulating 
layer  combinations.  Also,  they  are  representative  of  the 
order  in  which  one  would  use  the  algorithms  included  in 
the  appendix  in  order  to  analyze  a  particular  coating 
configuration.  That  is,  assuming  a  certain  material  is 
available  and  has  a  measured  electrical  conductivity, 
how  does  one  spatially  arrange  it  for  optimum 
effectiveness? 

First  of  all,  in  addition  to  selecting  a  coating 
material,  we  must  select  a  suitable  thermal  insulator. 

The  four  primary  areas  of  concern  in  selecting  an 
insulator  are  its  electrical  characteristics  which 
include  permittivity  and  electrical  conductivity  (or 
loss  factor),  its  thermal  conductivity,  and  its 
mechanical  applicability.  Extensive  experimentation 
was  conducted  using  Plexiglas,  styrene  plastic, 
styrofoam,  urethane  foam,  neoprene  rubber,  glass, 
phenolic,  paraffin  wax,  paper,  etc..  It  was  found  that 
the  foams  (styrofoam  and  urethane  foam)  exhibited  superior 
electrical  characteristics  over  the  .other  materials. 

Their  loss  factors  are  very  low  and  both  have  a 
relative  permittivity  very  near  1;  thus,  the  microwaves 
hardly  know  the  insulators  exist.  In  addition  they  are 
excellent  thermal  insulators.  Therefore,  where 
mechanical  application  is  relatively  straight  forward. 
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either  of  the  foams  will  perform  excellently.  Urethane 
foam  is  impervious  to  nearly  all  chemicals  and,  therefore, 
nearly  any  coating  may  be  used  without  reacting  with  it. 
Styrofoam,  on  the  other  hand,  is  a  styrene  and  exhibits 
disastrous  results  when  used  with  nearly  any  petroleum 
based  coating.  This  is  especially  true  with  toluene 
which  is  the  carrier  for  the  aquadaq  coating.  It  was 
found  that  styrofoam  performed  very  well  with  paraffin/ 
carbon  coatings  or  by  chemically  shielding  it  from  a 
petroleum  based  product.  The  shield  may  be  a  very  thin 
layer  of  paraffin  wax, stick-on  Mylar,  or  stick-on  paper. 
For  flat  surfaces,  there  is  a  %  inch  foam  sheet 
available  at  most  art  supply  stores  with  a  paper  coating 
on  each  side.  It  is  extremely  easy  to  work  with. 
Inexpensive,  and  comes  in  sheets  up  to  3x4  feet  square. 
When  working  on  complex  surfaces,  the  neoprene  rubber 
exhibited  the  most  flexibility  from  an  ;'pplications 
point  of  view;  however,  it  has  a  very  large  relative 
permittivity.  The  permittivity  problem  may  be  partially 
overcome  by  using  neoprene  foam  which  is  commercially 
available  in  thicknesses  from  1/16  to  1/2  inch.  Even 
though  the  permittivity  for  neoprene  foam  was  not 
measured  it  should  be  significantly  less  than  the  value 
of  31  measured  for  solid  neoprene.  This  is  the  same 
foam  used  in  diver's  wet  suits  and  should  also  be 
relatively  inexpensive.  It  is  easily  applied  to  almost 
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any  grease  free  surface  using  a  typical  rubber  contact 
cement.  It  is  also  impervious  to  nearly  all  chemicals. 

Having  chosen  a  coating  type  and  insulator,  the 
final  question  that  must  be  answered  is  related  to  the 
thickness  of  each.  Using  the  algorithm  included  in  the 
Appendix,  we  may  generate  a  set  of  three  contour  plots 
similar  to  the  ones  presented  earlier  in  this  chapter. 

The  first  would  be  only  with  the  coating  and  insulator 
and  no  substrate  present.  This  will  give  us  an  idea  of 
the  differential  temperatures  that  might  be  observed  on 
the  coating  if  there  were  zero  surface  currents  located 
on  the  substrate.  Next,  we  consider  the  coating/ 
insulator/substrate  combination.  The  temperatures 
illustrated  would  represent  the  AT^  that  might  be 
expected.  This  is  the  same  AT^  discussed  in  the  latter 
part  of  Chapter  VII.  In  this  case  the  surface  currents 
on  the  substrate  are  not  zero  but  something  between 
zero  and  The  third  step  is  to  numerically  subtract 

the  above  two  results  and  display  the  difference  as  a 
contour.  This  is  probably  the  most  important  step  since 
it  provides  the  insight  required  to  choose  the  optimum 
thickness  for  both  the  coating  and  insulator.  It  allows 
us  to  choose  thicknesses  that  will  provide  the  greatest 
temperature  difference  between  the  case  of  no  substrate 
(J*0)  and  infinite  substrate 
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Illustrative  of  the  above  process  are  two 
coating  examples  that  have  provided  excellent  results. 

The  first  is  three  layers  of  aquadaq  (315  mhos/m), 
yielding  a  total  thickness  of  approximately  15  jim, 
placed  on  a  5.1  mm  thick  sheet  of  paper  covered  foam. 

This  was  the  scheme  used  in  the  resolution  experiments 
in  Chapter  VII.  The  second  was  a  paraff in/carbon 
conductive  coating  used  by  Sega  approximately  1.5  mm 
thick  (=2.0  mhos/m)  and  placed  on  a  layer  of  styrofoam 
12,5  mm  in  thickness.  There  were  some  difficulties  with 
this  scheme  in  that  the  coating  thickness  tended  to  vary 
greatly  even  after  machining.  This  resulted  primarily 
from  foam/paraffin  expansion  problems;  therefore,  if  the 
paraff in/carbon  coating  is  to  be  used,  it  should  be  kept 
thin  (preferably  less  than  1  mm). 

In  summary  then,  for  the  case  of  a  coating/ 
insulator/substrate  experiment,  attempt  to  use  a  foam 
as  the  insulator:  preferably  styrofoam,  urethane,  or 
last  choice,  neoprene  foam.  Use  a  coating  that  is  easily 
applied  to  a  uniform  thickness.  In  this  case  differential 
surface  temperature  is  almost  directly  proportional  to 
coating  thickness  so  uniformity  is  paramount.  Typically, 
the  electrical  conductivity  will  be  a  fixed  value; 
however,  for  a  spray  any  value  from  10  to  500  mhos/m 
will  probably  work  well  (A  conductivity  between  20  to 
100  mhos/m  would  be  optimum.).  Execute  the  above 
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mentioned  computer  program  and  determine  the  optimum 
thicknesses  for  the  coating  and  insulator.  Apply  the 
coating/insulator  to  the  substrate  and  after  suitable 
curing  the  experiment  should  be  ready  for  thermographic 
observation  in  the  microwave  field. 

Most  of  the  preceding  discussion  has  centered 
predominantly  on  the  coating/insulator/substrate  problem; 
however,  it  should  be  pointed  out  again  that  this  is  not 
the  only  coating  technique.  If  we  are  interested  only  in 
the  geometrical  shape  of  a  particular  object,  rhe  coating 
problem  is  much  simplified.  The  idea,  as  was  mentioned 
earlier,  is  to  construct  a  model  of  foam  and  then  coat 
it  with  2  to  3  skin  depths  of  conductive  coating.  In 
this  case  it  can  be  shown  that  the  surface  temperature 
is  independent  of  coating  thickness  and  only  a  function 
of  the  surface  currents  present;  thus,  the  coating 
application  need  not  be  critical.  Clearly,  if  the 
coating  has  no  conductivity  (a  dielectric)  or  infinite 
conductivity  (perfect  conductor),  there  will  be  no 
energy  absorption  and,  thus,  no  surface  temperature 
increase.  The  primary  problem  then  is  to  determine  the 
optimum  conductivity.  Chapter  V  provided  the  answer 
for  us  by  referring  to  Figure  19.  There  it  can  be  seen 
for  an  input  power  of  10  mW/cm  and  at  a  frequency  of 
2.45  GHz  that  the  maximum  conductivity  we  may  have, 
and  still  have  a  suitable  temperature  increase  (1°  K) , 


137 


is  750  mhos/m.  If  we  need  to  consider  different  microwave 
parameters,  another  three  dimensional  plot  may  be  generated 
with  the  computer  program  included  in  the  Appendix.  In 
this  coating  scheme  we  must  keep  in  mind  that  we  want  the 
electrical  conductivity  as  large  as  possible  so  that  it 
will  more  closely  approximate  a  perfect  conductor.  If  we 
absolutely  need  a  conductivity  greater  than  750  mhos/m, 
we  must  increase  our  input  power  correspondingly  in  order 
to  observe  acceptable  heating. 

In  conclusion  we  may  reiterate  some  of  the 
important  aspects  of  this  paper.  Initially,  a  theoretical, 
one  dimensional  model  is  presented  which  couples  the  N- 
Layer  electromagnetic  problem  with  the  thermodynamic 
problem;  hence,  we  have  a  model  which  predicts  equilibriur 
surface  temperatures  resulting  from  electromagnetic 
absorption  in  a  system  of  N  layers.  This  is  provided 
as  an  analysis  tool  for  the  engineer  interested  in 
investigating  different  coating  schemes.  The  computer 
program  which  solves  this  problem  is  included  in  the 
Appendix.  This  model  is  subsequently  verified  experi¬ 
mentally  on  particular  multi-layered  models.  The 
applicability  of  the  small  sample  verification  process 
was  discussed  next.  It  was  shown  that  a  small  sample 
may  be  a  valid  approximation  for  the  one  dimensional 
model.  It  was  also  shown  that  nearest  neighbor 
considerations  are  minimized  for  coatings  that  are 
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placed  near  the  substrate  surface.  As  a  rough  rule  of 
thumb,  it  was  shown  that  the  coating  resolution  may  be 
approximated  by  the  value  a/ it,  where  "a"  is  the  insulator 
thickness.  It  was  also  shown  that  the  differential 
temperature  observed  on  the  coating  of  an  electrically 
large  object  may  be  thought  of  as  a  steady  state  term 
resulting  from  the  one  dimensional  solution  plus  a 
higher  order  term  resulting  from  the  finite  boundary 
conditions  present.  As  a  result  of  this  analysis,  the 
experimental  engineer  may  be  allowed  to  scale  his  infrared 
results  by  an  additive  constant  in  order  to  "calibrate” 
the  thermographic  system.  In  addition  to  the  above,  two 
particular  coating  schemes  are  illustrated  that  provided 
reasonable  results;  these  were  the  carbon/paraffin/ 
styrofoam  and  aquadaq/styrofoam  coatings.  Lastly, 
the  Appendix  includes  information  concerning  techniques 
for  measuring  electrical  permittivity  and  electrical 
conductivity  as  well  as  most  of  the  major  computer 
routines  used. 

Finally,  it  is  important  to  reiterate  that  this 
entire  development  was  classical  in  nature  and  thus  only 
allowed  for  joule  heating  in  the  coating.  Other 
absorption  mechanisms,  such  as  rotational  coupling,  may 
prove  equally  or  more  effectiv^e  in  coating  design  in  the 
future.  The  thermographic  detection  of  induced  surface 
currents  is  only  beginning  to  demonstrate  its  full 
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capability.  The  benefits  of  such  a  scheme  should  prove 
invaluable  to  the  systems  design  engineer  working  in  the 
ever  increasing  electronic  warfare  environment. 
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APPENDIX  A 

ELECTRICAL  CONDUCTIVITY  MEASUREMENTS 

The  purpose  of  the  electrical  conductivity 
measurements  was  to  determine  the  value  of  the  electri¬ 
cal  conductivity  of  the  various  coatings  used  in  the 
infrared  current  measurement  schemes . 

The  most  likely  source  of  inconsistency  in  the 
measurement  process  resulted  from  the  method  used  to 
attach  or  input  the  curreut  to  the  test  specimen.  The 
two  point  measurement  technique  used  in  this  testing 
is  recognized  by  the  American  Society  of  Testing  and 
Materials  (ASTM)  as  the  most  precise  method  for  deter¬ 
mining  conductivity.  However,  care  must  be  exercised 
to  assure  that  the  electrical  current  passed  through 
the  specimen  is  uniformly  distributed  over  the  entire 
cross  sectional  area  of  the  test  specimen. 

Tests  were  conducted  using  two  possible  specimen 
configurations.  The  first  consisted  of  a  10cm  x  10cm 
square  of  the  material  placed  on  a  plexiglas  substrate 
(Fig  96  )  . 
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Figure 
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Coating  Sample  Cross  Section 
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This  arrangement  was  used  predominately  for  thin  coat¬ 
ings  that  were  sprayed  on  the  substrate  from  aerosol 
cans;  thicknesses  varied  from  10-100  microns.  The 
other  configuration  which  was  used  with  "moldable" 
materials  such  as  carbon  impregnated  paraffin  consisted 
of  solid  disks  2.54  cm  in  diameter  and  approximately 
7mm  thick  (Fig  97) . 


Silver  Paint 


Silver  Paint 


Figure  97  :  Drawing  of  a  Typical  Disc  Sample 

In  both  cases  a  silver  metal  based  lacquer  was  used  to 
make  final  contact  with  the  specimen  thus  helping  in¬ 
sure  a  uniform  electrical  contact. 
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The  current  connections  for  the  lOctn  square 
were  attached  to  each  end  of  the  specimen  by  means  of 
a  copper  foil  strip  embedded  in  the  silver  lacquer. 
Connections  for  the  disk  were  made  by  sandwiching  it 
between  copper  plates  which  were  subsequently  clamped 
in  three  places  with  wave  guide  clamps.  See  Figures 
98  and  99. 


Figure  98:  Square  Sample  Electrical  Connections 


Figure  100:  Schematic  of  Conductivity  Measurements 
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This  test  arrangement  Is  the  ASTM  two-point  method. 

The  following  test  equipment  was  used  in  the  test 
arrangement : 

Nomenclature  Identification 

Power  Supply  Regulated  D.C.  Power  Supply 

0-50  VDC,  0-1. 5A 

Kepco  Mfg  Co,  Flushing,  NY 

Ammeter  Digital  Multimeter 

Hewlett-Packard  Model  3466A 

Voltmeter  Digital  Multimeter 

Hewlett-Packard  Model  3466A 

Micrometer  Metric  Micrometer 

Central  Scientific  Co, 
Chicago ,  IL 

The  conductivity  (o)  for  each  of  the  samples 
was  calculated  from  the  following  equation: 

0  »  il 

VA 

where  i  -  total  current  through  specimens 
1  «  length  of  current  travel 
V  -  voltage  across  sample 
A  -  area  through  which  current  < ravels 
Conductivity  versus  temperature  was  also 
investigated  for  the  10cm  x  10cm  samples  with  an 
aquadaq  coating.  An  AGA  Thermovision 680  infrared 
camera  was  used  to  measure  the  steady  state  surface 
temperature  of  the  sample  for  various  current  inputs. 
See  Figure  101  for  a  plot  of  o  versus  temperatures  (u) . 


ELECTRICfiL  CONDUCTIVITY,  MHOS' 
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COATING  TEMPERATUPE  in  DEGREES  t ELVIM 


Figure  101:  Electrical  Conductivity  Versus  Absolute 
Temperature  (Aquadaq  Coating) 
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Using  the  technique  outlined,  the  measured 
conductivity  for  a  composite  mixture  of  paraffin/ 
carbon  was  measured  for  various  mixing  ratios.  These 
plots  simplify  the  process  of  designing  a  particular 
conductivity  material  since  the  particular  mixing 
ratios  are  illustrated.  See  Figures  102  and  103  for 
a  linear  and  log  plot  of  sigma.  The  solid  line  is 
a  plot  of  an  empirical  model  of  sigma  as  given  in 
Figure  103. 


Sigma  in  Mhos 
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EXPEPIMENTRL  ELE.CTPICRL  CONDUCT  IV I  T  i 
^.00E+01 


• - . 5*0' . '6fe"”"7fe . '4\h' . 'sfe . i'^0 

Conductor  Weight  by  Per  Cent 


Figure  102:  Electrical  Conductivity  Versus  Carbon/Paraffin 


Mixing  Ratios 
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EXPEPlriENTRL  EL.  ECFRICRL  CONDUCTlVni 


Ha»irnum  b  i  gm  a-*!  .  6.i£! +0 1 
P-+6.845E+1  1 
O=+3.468E^0^ 


Conductor  Weight  By  Per  Cent 


Figure  103:  Electrical  Conductivity  Versus  Carbon/Paraffin 
Mixing  Ratios  with  Empirical  Model 
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APPENDIX  B 

ELECTRICAL  PERMITTIVITY 

Electrical  permittivity,  e,  is  an  exceedingly 
important  variable  in  the  analysis  of  effects  of 
microwave  absorption  in  a  multi-layered  system.  It  is 
typically  frequency  dependent  and  therefore  the  value 
used  must  be  measured  for  each  different  operating 
condition.  Fortunately,  it  does  not  vary  greatly;  the 
table  below  lists  values  taken  from  the  Chemical 
Rubber  Company's  Handbook  of  Chemistry  and  Physics, 
60th  Edition,  for  some  common  materials. 


Table  1:  Relative  Permittivity  for  Various  Dielectrics 


Sodium 

Light 

iMHz 

100  Mhz 

Specific 

Heat 

Paraffin 

2.0 

2-2.5 

* 

* 

Glass 

2. 3-3. 6 

4.0 

* 

* 

Rock  Salt 

2.3 

* 

* 

* 

Gelatin 

2.3 

* 

★ 

* 

Quartz 

2.3 

* 

* 

* 

Nylon  (66) 

2.3 

3.3 

3.2 

.4 

Polyethylene 

2.3 

2.3 

2.3 

.55 

Methylmethacrylate 

(Plexiglas) 

2.2 

2.8 

* 

.35 

Polystyrene 

2.6 

2. 5-2. 7 

2.6 

.32 

Silicon  Rubber 

* 

3. 1-3. 2 

* 

* 

Porcelain 

* 

6-8 

* 

* 

♦Value  not  available 
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Our  permittivities  were  measured  in  the  labora¬ 
tory  at  a  frequency  of  10  GHz  for  styrofoam,  window  glass, 
paraffin,  styrene,  phenolic,  plexiglas,_  and  neoprene. 
Figure  104  illustrates  the  experimental  arrangement. 

Test  equipment  used  in  this  measurement  scheme  included 
the  following; 

Nomenclature  Identification 

10  GHz  Source  ED- SET  MARK  2 

Sargeant-Welch  Scientific  Co. 
Skokie,  IL 

10  GHz  Detector  " 

Microwave  Bench  ” 


Figure  104  :  Schematic  Arrangement  of  Permittivity 


Measurements 
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The  procedure  was  to  measure  the  reflectance 
for  a  given  sample  at  various  angles  of  incidence.  E 
was  polarized  perpendicular  to  the  plane  of  incidence. 
Because  of  limited  sample  sizes,  incident  angles  were 
limited  to  between  20  and  60  degrees  with  measurements 
taken  in  5  degree  increments.  The  measured  value 
versus  angle  was  then  plotted  on  a  graph  on  which  theo¬ 
retical  curves  of  reflectivity  versus  permittivity 
had  been  drawn  for  a  particular  sample  thickness  and 
incident  angle.  The  permittivity  is  where  the  measured 
reflectance  value  intersects  the  particular  curve. 
Figures  105  to  ill  illustrate  the  plots  for  styrofoam, 
glass,  paraffin,  styrene,  phenolic,  Plexiglas,  and 
neoprene  respectively.  (Only  10  degree  increments 
are  plotted  for  illustration  purposes.) 


Reflectivity  in  Per  Cent 
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l-LRYER  ELECTROMAGNETIC  WAVE  INTERACTION 


Figure  105:  Permittivity  Plot  for  Styrofoam  at  10  GHz 


Pe -f  1  ec  t  1  V  1 1 >  >  ri  Per  C en 
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Figure  106 


Permittivity  Plot  for  Window  Glass  at  10  GHz 


Re-f  lect  1  vity  in  Per  Cent 
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1 -LAYER  ELECTROMAGNETIC  NAVE  INI  ERR'  I  L'N 


Figure  107  :  Permittivity  Plot  for  Paraffin  at  10  GHz 


Reflectivity  in  Per  Cent 
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Samp  I e  Cpc  <  ) on  ( r  e ) at i ve  ^ 


l -LAYER  ELECTROMAGNETIC  WAVE  INTERACTION 


Figure  108:  Permittivity  Plot  for  Styrene  at  10  GHz 


Reflectivity  in  Per  Cent 
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S«mp  I  •  ^  on  (  r  e  1  at  i  va  > 


l -LAYER  ELECTROMAGNETIC  WAVE  INTERACTION 


Figure  109  ;  Permittivity  Plot  for  Phenolic  at  10  GHz 


Re-f  lect  I  V  I  ty  Per  Cent 


S  1  «  Ep  &  I  1  on  '  r  «  I  i  VC  ) 

l-LRYER  ELECTROnflGNETIC  WAVE  INTERACTION 


Figure  110  :  Permittivity  Plot  for  Plexiglas  at  10  GHz 


a  9.0  t7'.0  25.0  33.0  4  1.0 

Sftmp I e  £p& I  1  on  ( ( e 1 i v«  ) 


2 --LAYER  ELECTROriAGrJETIC  WAVE  INTERACTION 


111  :  Permittivity  Plot  for  Neoprene  at  10  GHz 


No  attempt  was  made  to  measure  permittivities 
for  any  materials  other  than  good  dielectrics.  Since 


Ohm's  law  was  included  explicitly  in  the  solution  of 
Maxwell’s  equations,  it  is  not  necessary  to  assume  a 
complex  permittivity  for  conductive  materials;  thus, 
the  free  space  value  is  acceptable. 

Table  2  lists  the  measured  values  of  rela¬ 
tive  permittivity  for  the  materials  considered.  The 
averaging  procedure  consisted  of  calculating  an 
average  value  and  a  standard  deviation  for  the  nine 
incident  angles.  Then  any  measured  permittivity  out¬ 
side  the  range  *  one  standard  deviation  from  the  average 
was  rejected  and  a  new  average  and  standard  deviation 
calculated.  These  are  the  values  that  appear  in  Table 
2  below. 

Table  2:  Measured  Relative  Permittivity  at  10  GHz 

Material  Average  Relative  Permittivity  (10  GHz) 


Styrofoam 

1.11 

(One 

data  point) 

Window  Glass 

14.2 

(St. 

Dev. 

s 

1.3) 

Paraffin 

2.55 

(St. 

Dev. 

s 

.14) 

Styrene 

4.28 

(St. 

Dev. 

= 

.08) 

Phenolic 

5.32 

(St. 

Dev. 

s 

.33) 

Plexiglas 

4.07 

(St. 

Dev. 

s 

.24) 

Neoprene 

31.1 

(St. 

Dev. 

= 

.83) 
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APPENDIX  C 

COMPUTER  PROGRAM 


The  name  of  this  program  is  "One-B".  It  is 
written  for  use  on  a  Hewlett-Packard  9845B  minicomputer. 
It  calculates  reflectivity,  transmissivity,  and  absorp¬ 
tivity  for  a  single  interface.  It  will  also  calculate 
infrared  emissivity  for  5.3  microns  wavelength.  Input 
variables  are  medium  two  electrical  conductivity,  index 
of  refraction,  and  magnetic  permeability  (relative)  as 
well  as  the  incident  frequency  in  gigahertz. 


!  THE  NAME  OF  THIS  PROGRAM  IS  ONE-B ! 

INTEGER  I, J,0,N,K,L,M 

DIM  Para<4,5>,Perp<:4,5>,  Theda<90>,Refpar<90),Refperf90>,  Jeffpa<90> 
DIM  Jeffpe<90) 

GCLEAR 
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The  name  of  this  program  is  "Basic".  It 
calculates  the  reflectivity,  transmissivity,  and 
absorptivity  for  an  N-layer  electromagnetic  interaction 
problem.  The  input  parameters  are  the  number  of  layers 
(not  to  exceed  10),  the  incident  frequency  in  gigahertz, 
the  incident  power  in  mW/cm  ,  and  the  individual  layer 
conductivities,  permittivities,  permeability,  and 
thickness . 


THE  NfiME  OF  THIS  PROGRfiM  IS  'Basic'  ! 

IT  CnLCULRTES  THE  REFLECTANCE,  TRANSMITTANCE,  ANB  ABSORBANCE  FOR  A  N-LAY 
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APPENDIX  E 
COMPUTER  program 

The  name  of  this  program  is  "SIGNIF".  It 
calculates  the  ratio  of  heat  lost  through  infrared 
radiation  to  that  lost  by  way  of  heat  convection  on 
a  vertical  flat  surface.  The  input  variables  are 
air  temperature  (C) ,  the  convection  exponent,  and 
plate  height  (m) . 
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APPENDIX  F 

COMPUTER  PROGRAM 

The  name  of  this  program  is  "Uo-3D".  It 
calculates  the  equilibrium  surface  temperature  for  a 
system  of  N  layers  in  the  presence  of  electromagnetic 
radiation.  It  plots  the  temperature  as  a  function  of 
the  first  layer  thickness  and  conductivity  on  a  three 
dimensional  plot.  The  input  variables  are  the  number 
of  layers,  the  incident  power,  the  incident  frequency, 
the  ambient  air  temperature,  the  convection  exponent, 
the  surface  emissivity,  the  vertical  flat  plate  height, 
the  incident  angle,  and  the  smallest  electrical 
conductivity  for  the  first  layer.  Additionally,  the 
inputs  for  the  various  other  layers  are  the  permittivity, 
conductivity,  permeability,  and  thickness. 
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APPENDIX  G 
COMPUTER  PROGRAM 

The  name  of  this  program  is  "Uthick".  It 
calculates  the  surface  temperature  of  a  system  of  N 
layers  as  a  function  of  the  thickness  of  any  two  of 
the  layers.  The  result  is  plotted  on  contour  plot 
of  constant  temperatures.  The  input  variables  are  the 
number  of  layers,  which  layers  to  vary  the  thickness 
of,  the  incident  power,  the  incident  frequency,  the 
incident  angle,  the  ambient  air  temperature,  the  plate 
height,  surface  emissivity,  the  convective  exponent, 
and  the  layer  characteristics  of  conductivity, 
permittivity,  permeability,  and  thickness. 


3ie  RSSIGN  «1  TO  "UoBIG 

320  GOTO  360 
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APPENDIX  H 
COMPUTER  PROGRAM 

The  name  of  this  program  is  "CONTOR".  It 
will  plot  and  smooth  contour  lines  for  an  arbitrary 
matrix  of  values.  It  will  do  linear  smoothing,  circular 
or  parabolic  smoothing.  The  matrix  data  is  read  from 
a  data  file  named  "X’Data".  Plot  labeling  is  accomplished 
manually  when  the  flashing  cursor  appears  on  the  9845 
CRT.  The  only  inputs  are  the  number  of  rows  and 
columns  of  the  input  matrix. 
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APPENDIX  I 

COMPUTER  PROGRAM 

The  name  of  this  program  is  "EPSLON".  It 
calculates  the  per  cent  of  reflected  to  incident  power 
and  plots  it  versus  relative  permittivity  for  the  layer. 
The  input  variables  are  the  number  of  layers,  the  maximum 
permittivity  to  be  considered,  the  incident  frequency, 
the  incident  angle,  and  the  measured  reflectivity. 


THE  NfiME  OF  THIS  PROGRAM  IS  'EPSLON 
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